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ABSTRACT 

“Green Bean” Galaxies (GBs) are the most [O III]-luminous type-2 active galactic 
nuclei (AGN) at 2 : ^ 0.3. However, their infrared luminosities reveal AGN in very 
low activity states, indicating that their gas reservoirs must be ionized by photons 
from a recent high activity episode - we are observing quasar ionization echoes. We 
use integral field spectroscopy from the Gemini Multi-Object Spectrograph to analyse 
the 3D kinematics, ionization state, temperature and density of ionized gas in the 
GB SDSS J224024.1-092748. We model the emission line spectrum of each spaxel as 
a superposition of up to three Gaussian components and analyse the physical prop¬ 
erties of each component individually. Two narrow components, tracing the velocity 
fields of the disc and an ionized gas cloud, are superimposed over the majority of the 
galaxy. Fast shocks produce hot (Tg ^ 20,000 K), dense (rig ^ 100cm“^), turbulent 
(<j ^ 600kms“^), [0111]-bright regions with enhanced [NII]/H(a and [SII]/H(a ratios. 
The most prominent such spot is consistent with a radio jet shock-heating the inter¬ 
stellar medium. However, the AGN is still responsible for > 82 per cent of the galaxy’s 
total [O III] luminosity, strengthening the case for previous quasar activity. The ion¬ 
ized gas cloud has a strong kinematic link to the central AGN and is co-rotating with 
the main body of the galaxy, suggesting that it may be the remnant of a quasar-driven 
outflow. Our analysis of J224024.1-092748 indicates that GBs provide a unique fossil 
record of the transformation from the most luminous quasars to weak AGN. 

Key words: galaxies: active - galaxies: evolution - galaxies: individual: SDSS 
J224024.1-092748. 


1 INTRODUCTION 

Extended ionized gas reservoirs around active galactic nuclei 
(AGN) are valuable probes of the mechanisms driving rapid 
changes in AGN accretion rates. Kiloparsec-scale regions of 
highly ionized gas (narrow line regions; NLRs) are excited 
by energetic outflows, driven by either radiation pressure (in 
luminous AGNs with high accretion rates; see e.g. |Alexan-] 
der et aL]|2010 ) or radio cores/jets (in radio-loud systems; 


Canalizo 2006 for a review). For this reason, EELRs have 


Rosario et al.|20i0 Mullaney et al.|2013| ) . The most 


luminous AGN can ionize gas to distances of tens of kilo- 
parsecs, forming extended emission line regions (EELRs). 
EELRs have significantly lower surface brightnesses than 
typical AGN NLRs, but require substantially more energy 


to ionize due to their large volumes (see Stockton, Fu & 
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been observed almost exclusively around radio-loud quasars 
dFu fc Stockton||2009 ). 

The small subset of EELRs found around galaxies with 
weak AGN are ionization echoes and provide fossil records of 
the rapid shutdown of luminous quasars. The long recombi¬ 
nation timescale in low density EELR gas allows it to remain 
highly ionized for tens of thousands of years after the quasar 
phase of the central engine has ceased. Hanny’s Voorwerp, 
a large and highly ionized gas cloud associated with the 
galaxy IC 2497 (z 0.05), was discovered serendipitously in 
the Galaxy Zoo survey ( [Lintott et al.|200^ Schawinski et al. 
2010 Keel et al.|2012b| ). The lack of current AGN activity in 
IC 2497 indicates that the luminosity of the AGN must have 
decreased by at least 2-3 orders of magnitude in less than 
the light crossing tim e from the nucleus to Hanny’s Voorw- 


erp (10^ years, see Lintott et al. 2009 Schawinski et al 


2010| [Keel et ^|2012b ). A targeted Galaxy Zoo follow-up 
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search revealed 8 further AGN ionization echoes at z < 0.1 


with similar properties to Hanny’s Voorwerp (Keel et ah 


2014). These objects provide unique insights into both the 
life cycle of AGNs and the physical processes driving rapid 
changes in AGN accretion rate. 

The descendants of the most luminous quasars are very 
rare and can therefore only be found by increasing the 
search volume to higher redshifts. The prototypical quasar 
ionization echo SDSS J224024.1-092748 (J2240 hereafter), 
at a redshift of z = 0.326, was discovered serendipitously 
in a Ganada France Hawaii Telescope (GFHT) image by 
Schirmer et al. ( 2013| ). A follow-up search using the Sloan 
Digital Sky Survey (SDSS) Data Release 8 archive pro¬ 
duced a sample of 29 ‘Green Bean’ galaxies (GBs) across 
two redshift ranges (0.19 < z < 0.35 and 0.39 < z < 0.69). 
The GBs are defined by their extreme [O III] A5007 fluxes 
(r ~ 18 mag at z ~ 0.25) and large angular extents (several 
arcseconds; Schirmer et al.|[^013 ). GBs occupy a volume in 
color-magnitude-size space that is populated by about 95 
per cent spurious detections which must be removed by vi¬ 
sual inspection. 

GBs are the most luminous AGN at their cosmic time. 
Their EELRs span 20 — 40 kpc in diameter and have typical 
[GUI] luminosities of \og{L[oiii] / erg s“^) = 42.7 - 43.6, 
making them ~ 10 times more luminous than the 1300 refer¬ 
ence Type-2 quasars at similar redshifts selected fromjRey^ 


et al. (2008) and Mullaney et al. (2013). Only half of the 


known GBs have radio luminosities larger than 10^"^ WHz“ 
compared to 98 per cent of all type-2 AGN at similar 
redshifts (Mullaney et al. 2013), demonstrating that their 
EELRs must be excited largely by radiation pressure. The 
EELRs of GBs lie well above the empirical size-luminosity 
relation (see e.g. Liu et al. 2013), indicating that their 


extreme luminosities may be a by-product of larger-than- 
average gas reservoirs (see e.g. Hainline et al.||2013 ). 

The most unusual property of GBs is that their 22/im 
WISE IR luminosities (L 22 nm) are 5 — 50 times lower than 
expected from L[o m] (| Schirmer et al.|| 


and L\o m] probe AGN activity ( jAsmus et al.[2011 


Ichikawa 

et al.|201^ . The infrared (IR) flux in active galaxies is emit- 


Both L 


22fj,m 


ted by heated dust within 10—100 pc of the SMBH, whereas 
the [O III] emission originates from the EELR on kilopar- 
sec scales. Therefore, the discrepancy between L 22 ^im and 
^[o III] iriay be best explained by a recent AGN shutdown, 
which due to the long gas recombination timescale has not 
yet propagated through the EELR. Another explanation is 
that the dusty torus emits disproportionally less IR light at 
a given AGN bolometric luminosity, due to e.g. structural 
differences in the torus (Hai Eu, private communication). 
However, it is not conceivable that the SDSS color selection 


adopted by Schirmer et al. (2013) could bias the luminosity 
of our sample enough to produce such an effect. Processes in¬ 
volving the orbital decay of binary AGN or the coalescence 
of two SMBHs cannot be responsible for the discrepancy, 
as the X-ray emissions involved are too low and short-lived 


( Hopkins Quataert|201Q) Tanaka, Haiman &: Menou|201Q 


khan et al.||2012 ). 

Preliminary analysis of new Ghandra data for 9 GBs 
reveals X-ray counts much lower than expected from their 
mid-IR luminosities in the limiting case of Gompton-thick 
(log[A^H] > 24.2 cm“^) absorbers. Either these AGN are 
deeply buried, or we are witnessing the descendants of lumi- 


Channel 

Blue 

Green 

Red 

Instrument 

GMOS-N 

GMOS-S 

GMOS-N 

Grating/Filter 

B600 + g 

B600 + r 

R831 -j- z -j- CaT 

Date obtained 

2012-10-08 

2012-08-27 

2012-10-18 

Xobs (nm) 

429 - 542 

559 - 688 

847 - 905 

Xrest (nm) 

323 - 409 

421 - 520 

638 - 683 

R 

2700 

2700 

7030 

Exposure time (s) 

4 X 1800 

6 X 1500 

4 X 1600 

Mean airmass 

1.17 

1.10 

1.19 


Table 1. Characteristics of CMOS integral field unit (IFU) ob¬ 
servations. 


nous quasars 10^“^ years (i.e. much less than the light cross¬ 
ing time of the EELR) after the shutdown of their quasar 
phases. The large angular extent of the EELRs allow them 
to retain a memory of the past X-ray histories of the AGNs, 
spanning up to ~ 10^ years. Our Ghandra observations will 
be described in detail in a future paper. 

Using photoionization models and resolved 3D optical 
spectroscopy, it is possible to extrapolate the incident X-ray 
flux at a given distance from the nucleus to infer the nu¬ 
clear X-ray luminosity for different times in the past. The 
EELRs of GBs are sufficiently bright and luminous that light 
curve reconstruction is possible. However, such reconstruc¬ 
tion requires a robust estimate of the contribution of shock 
ionization, and a detailed understanding of the formation, 
composition and physical state of the gas components in 
the EELR. In this paper we present deep optical integral 
held spectroscopy of the GB prototype J2240, spanning the 
330nm—680nm optical rest-frame wavelength range. We ex¬ 
plain our observations, data reduction procedure and data 
processing techniques in Section and decompose the gas 
into different kinematic components in Sectionj^ We present 
a detailed 3D spectroscopic analysis of the velocity held, di¬ 
agnostic line ratios, ionization parameter, and electron tem¬ 
perature and density in Section ^ and discuss the origin 
of the EELR in Section |5l We summarise our results and 
present our conclusions in Section 


2 OBSERVATION AND DATA PROCESSING 
2.1 Gemini observations 


We obtained deep Gemini Multi-Object Spectrograph 
(GMOS) observations of J2240 under observing programmes 
GS-2012B-Q-26 and GN-2012B-Q-226 (PL Schirmer). The 
two GMOS instruments (GMOS-S and GMOS-N) have 
four observing modes: imaging, long-slit spectroscopy, 
multi-object spectroscopy and integral held spectroscopy 


(Allington-Smith et al. 


2002 


Hook et al. 2004). Our in¬ 


tegral held unit (lEU) observations cover three spectral 
ranges: blue (B600 + g), Xobs = 429 - 542 nm, R 2700), 
green (B600 + r, Xobs — 559 - 688 nm, R 2700) and red 
(R831 + z + GaT, Xobs = 847 - 905 nm, R 7030). The 
wavelength coverage, spectral resolution, exposure times and 
dates and mean airmasses of our observations are sum¬ 
marised in Table [T] The observations were carried out in 
IQ 70 conditions, corresponding to a Point Spread Eunction 
(PSE) with a Eull Width at Half Maximum (EWHM) of 
approximately 0.5 - 0.6 arcsec (5-6 pixels) across all three 
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spectral channels. The angular extent of the target over 
which we conduct spectroscopic analysis is 2.5" x 5", within 
which there are 32-50 separate resolution elements. 

We observed spectrophotometric standard stars for flux 
and wavelength calibration of all three spectral channels: 
BD +28 4211 on 2012-10-05, LTT 9239 on 2012-10-31 and 
EG 131 on 2012-08-29, and Wolf 1346 on 2014-04-03 for the 
blue, green and red channels respectively. We also obtained 
the other standard CMOS baseline calibrations including 
night-time flats, twilight spectra, day-time arcs & biases. 


2.2 Data reduction 


Our data were reduced with the UrekeQl.O version of IRAF 


2.16 (Tody 1986 1993), using our own modified version 


of Gemini’s GMOS IRAF package, which we refer to as 
ifudrgmos, along with our PyFU datacube manipulation 
scripts. Our modifications are described in the Appendix. 
Both our ifudrgmos and PyFU packages are available to 
other investigators through the Gemini Data Reduction User 
Forunjf] For our GS-2012B-Q-26 and GN-2012B-Q-226 data 
we used revisions 79 and 114, respectively, of the Subversion 
repository containing these packages, which are functionally 
the same as the versions on the forum (since the numbering 
tracks changes to multiple packages). 

Our top-level processing sequence - based on the stan¬ 
dard GMOS data reduction - was as follows: bias and over¬ 
scan subtraction; creation and association of bad column 
masks; interpolation over bad columns using continuum fits; 
tracing of fibre spectra and gap regions using flat field ex¬ 
posures from the calibration unit, GGal; modelling and sub¬ 
traction of scattered light in scientific and flat exposures 
based on the count levels in gaps between blocks of fibres; 
identification of cosmic rays; cleaning of cosmic rays using 
interpolated replacement values; mosaicking of the GMOS 
detectors into a single array; extraction of fibre spectra 
to one dimension; wavelength calibration; normalization of 
GGal flats; division of scientific and standard star data by 
the flat field spectrum; rectification to linear wavelength co¬ 
ordinates; sky subtraction using the background field (sep¬ 
arately for each IFU slit); integration over the IFU field 
to obtain standard-star spectra; flux calibration with refer¬ 


ence to IRAF’s tabulated data from Bessell (1999), Hamuy 


et al. (1994) and Massey et al. (1988); correction for sub- 
A spectral flexure using sky lines; resampling the spectra 
to 3-dimensional ‘data cubes’; spatial registration (by cross¬ 
correlation) and co-addition of the data cubes for individual 
exposures within each wavelength channel; and re-binning 
of the co-added data cubes in logarithmic wavelength incre¬ 
ments, resulting in constant radial velocity increments (see 
Sect. 2.2.8). The last few processing steps were also repeated 
identically without sky subtraction, providing a measure of 
the final spectral resolution. 

Variance and data quality arrays were propagated 
throughout the process, overhauling a number of limitations 
and inaccuracies in the existing GMOS code. As usual, ran¬ 
dom measurement errors were initially estimated from the 
bias-subtracted data values in electrons (assuming a Poisson 


^ http://ssb.stsci.edu/ureka 
^ http://drforum.gemini.edu 


distribution) and the nominal detector read noise. The cal¬ 
culations do not include systematics, nor contributions from 
the high-signal-to-noise flat field and standard exposures. At 
each step, the variance estimate from the previous one was 
propagated using the standard uncorrelated error calcula¬ 
tions corresponding to the arithmetic operations performed 
on the main data array. When resampling, the variance ar¬ 
rays were interpolated in the same way as the scientific data, 
ignoring any minor smoothing and covariance introduced on 
small scales (which are zero in the ideal limit of sine inter¬ 
polation) . The greater covariance introduced by a change in 
pixel scale when creating data cubes was dealt with over¬ 
all using a simple scaling factor (see Section A8), to avoid 
having to track redundant, higher-dimensional information. 
Separate data quality bit planes were kept for detector de¬ 
fects, chip gaps and cosmic rays, allowing each to be treated 
at the most appropriate processing stage and in the corre¬ 
sponding dimensions. 

Differences in our data reduction from the standard 
GMOS procedure are described in the Appendix. Follow¬ 
ing publication, we will also provide an example GL script 
for this process (and the associated software improvements) 
through the Gemini forum. 


2.3 Data cube alignment and continuum fitting 


We fit the stellar continuum and nebular emission spectra 
of each spaxel across all three spectral channels simultane¬ 


ously (see Section 2.4), and therefore it is critical that the 
relative spatial alignment between the data cubes is accu¬ 
rate to within one pixel. We spatially align the data cubes 
by comparing flux maps for the central 75 kms“^ of the 
[OII] AA3727, 3729, [O III] A5007 and Ra A6563 emission 
line profiles (see Figure [^. Gontours trace the flux level in 
the Ha slice. The original data cubes (top row) have a rel¬ 
ative offset of a few pixels in the ^/-direction. However, this 
was accurately corrected for by eye using the spatial struc¬ 
ture in the flux maps. 

The stellar continuum is modelled as a linear combi¬ 
nation of stellar templates using the penalized pixel fitting 
routine (PPXF; [Cappellari Emsellem 2004). We adopt the 
high spectral resolution (AA — 0.3 A) stellar evolutionary 
synthesis models of Gonzalez Delgado et al. (2005), using 
isochrones from Bertelli et al.| ( |1994| ); [Girardi et al. ( |2000| ), 
and Girardi et al. ( |2002| . The templates cover 3 metallicities 
(Z = 0.019, 0.008, 0.004) and 74 ages (4 Myr < t < 17 Gyr) 
with a 0.06 log time sampling. The principal aim of the stel¬ 
lar population fitting is to correct for Balmer absorption 
produced by old stellar populations, and therefore the exact 
choice of stellar templates will have a minimal impact on 
our results. We subtract the fitted stellar continuum cubes 
from the full data cubes to produce emission line cubes. 


2.4 Emission line fitting 

3D analysis of emission line fluxes and emission line ratios 
is one of the most powerful probes of the ionizing sources 
and kinematics of emission line galaxies. Emission line lu¬ 
minosities probe the total energy output of individual ioniz¬ 
ing sources, whilst emission line ratios provide information 
regarding the temperature, density, metal abundance and 
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[Oil] 


[O III] 


Ha 



Figure 1. Top: Flux maps for the central 75 kms“^ of the [Oil] AA3727, 3729, [O III] A5007 and Ha A6563 emission line profiles in the 
original data cubes. Contours trace the flux level in the Ha slice. There is a relative offset of a few pixels in the y-direction between the 
cubes. Bottom: The same flux maps extracted from the aligned data cubes. The spatial variation in the flux levels is very similar across 
all slices; indicating that the cubes are well aligned. 


excitation state of the line-emitting gas. The wavelength 
centroid and shape of the emission line profiles also encode 
detailed information about 3D kinematic structure. 

The emission line spectra of galaxies can be analysed 
either parametrically or non-parametrically, with advan¬ 
tages and disadvantages associated with both methods. Non- 
parametric emission line flux measurements are made by 
integrating the flux within specific wavelength bandpasses. 
The gas velocity is estimated by comparing the intrinsic 
wavelength of an emission line to the wavelength where the 
flux peaks, and the velocity dispersion is estimated using 
the FWHM of the emission line. The principal advantage of 
non-parametric analysis is that no assumption is made re¬ 
garding the shapes of the line profiles. This allows for accu¬ 
rate measurement of the total flux in isolated emission lines; 
however, blended or doublet emission lines (e.g. [NII] A6548 
+ Ha A6563 + [Nil] A6583, [ S II] AA6717, 6731) c omplicate 
the analysis significantly (see [Schirmer et al.|2Q13| for an ex¬ 
ample of deblending these lines). 

Parametric analysis relies on the assumption that the 
observed emission line profiles can be reproduced using a lin¬ 
ear combination of analytic line profiles (most often Gaus¬ 
sian). Any number of Gaussian components can be fit to 
each emission line, and each component is parametrized by 
its velocity (relative to the systemic velocity of the galaxy), 
velocity dispersion and peak flux. Parametric emission line 
fitting codes overcome difficulties associated with blended 
lines and doublets by fitting the profiles of all emission lines 
simultaneously. Parametric analysis is most suited to sys¬ 
tems where multiple kinematic components are superim¬ 
posed along the line of sight and all contribute significantly 


to the total line emission of the galaxy. A fitted spectrum is 
produced for each kinematic component, allowing the phys¬ 
ical properties of different line-emitting regions to be anal¬ 
ysed independently. 

The kinematic complexity of J2240 (revealed by 
VLT/XSHOOTER long-slit data; see Schirmer et al.||2013 ) 
suggests that parametric emission line analysis would pro¬ 
vide significant insight into the ionizing sources of this 
system. We use a modified version of the IDL emission¬ 
line fitting toolkit LZIFU (Ho et ah, in prep, see Section 
4.1 of Ho et al. 2014 for a brief description) to measure 
the fluxes, velocities and velocity dispersions for a suite of 
14 temperature-, density- and excitation-sensitive emission 
lines ([OH] AA3726, 3729, [Ne HI] A3869, Hy, [O HI] A4363, 
Hell A4686, H/3, [OIII] A4959, [OIII] A5007, [Nil] A6548, 
Ha, [Nil] A6583, [S H] A6717 and [S H] A6731). Unfortu¬ 
nately, the [O I] A6300 line lies at a redshifted wavelength of 
8354 A and therefore does not fall within any of our selected 
spectral channels. LZIFU uses the Levenberg-Marquardt 
least-squares fitting technique (MPFIT; [Markwardt 2009) 
to model emission line profiles as a superposition of up to 
three Gaussian components. The reduced-x^ values of the 1, 
2 and 3 component fits are compared to determine the num¬ 
ber of statistically significant components in each spaxel. We 
require the Ha and [O III] fluxes of each significant compo¬ 
nent to be detected at the Scr level. All emission line maps 
shown in this paper only include spaxels detected at the So¬ 
le vel. 


The flexibility of LZIFU allows it to effectively reproduce 
a large variety of complex line profiles across J2240, as illus¬ 
trated in Figure This schematic shows the total reduced 



















J224024-1-092748: Quasar ionization echo 5 



Figure 2. Schematic showing reduced spectra (black) and emission line fits (individual kinematic components shown in red, blue and 
green) as well as the total fits (brown) for eight different regions in J2240. Red numbers in the top right hand corner of each sub-plot 
correspond to the numbered locations in the central image of the galaxy. The two upper panels in each sub-plot show spectra zoomed 
in on the + [OIH] complex (4800-5100 A; left hand side) and the [Nil] + Ha and [S II] complexes (6450-6800 A; right hand side). 
The two lower panels show the fitting residuals for each spectral region. The fluxes are normalized to the peaks of the [O III] and Ha 
lines in each spectral range individually. The fitting residuals are on the order of a few percent, indicating that our parametric emission 
line fitting successfully reproduces the observed spectra in J2240. 
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spectra (black) and emission line fits (individual kinematic 
components shown in red, blue and green) as well as the 
total fits (brown) for eight different regions in J2240. Red 
numbers in the top right hand corner of each sub-plot corre¬ 
spond to the numbered locations in the central image of the 
galaxy. The two upper panels in each sub-plot show spectra 
zoomed in on the H/3 + [OIII] complex (4800-5100 A; left 
hand side) and the [NII] + Ha and [S II] complexes (6450- 
6800 A; right hand side). These spectral regions contain the 
strongest emission lines in the wavelength range probed by 
our data and therefore demonstrate the complexity and vari¬ 
ety of the line profiles most effectively. The two lower panels 
show the fitting residuals for each spectral region. The fluxes 
are normalized to the peaks of the [O III] and Ha lines in 
each spectral range individually. The central panel shows the 
intensity within the central velocity slice of the [O III] A5007 
emission line. 


In referring to different regions of the galaxy, we adopt 
the nomenclature of Schirmer et al. (2013) whereby the 
south western region of the system is labelled the ‘cloud’ 
and the north eastern side the ‘disc’. We also note the pres¬ 
ence of two prominent emission peaks in the cloud region of 
the galaxy, which we label the northern and southern [O III] 
peaks. 


We observe a range of line profiles across the system. 
We note that the red channel is much noisier than the 
green channel due to the prominence of skyline residuals at 
large observed wavelengths (e.g. Xna^obs ~ 8700A). Near the 
outer edge of the cloud (region 3), the emission lines are very 
narrow and only a single kinematic component is present. 
The [OIII], Ha and [Nil] A6583 emission lines have large 
EWs (247, 47 and 18 A respectively), but [S H], [Nil] A6548 
and H/3 are barely detected above the noise. The remainder 
of the probed regions lie closer to the main body of the 
galaxy and therefore the [Nil], Ha, [OIII], H/3 and [S H] 
lines all have large EWs (ranging from tens to hundreds 
of Angstroms). The emission line fitting indicates that the 
spectra of these regions are best fit by a superposition of 
two narrow components with a small (~ 50 kms“^) veloc¬ 
ity separation, plus a prominent broader component whose 
velocity relative to the narrow components varies between 
different regions of the galaxy. 


The contribution of the stellar continuum emission to 
the total observed flux is negligible across the majority of 
the system (< 3 per cent of the peak Ha flux and <0.7 per 
cent of the peak [OIII] flux). We detect significant stellar 
continuum emission (7.7 per cent of the peak Ha flux and 
1.9 per cent of the peak [OIII] flux) in region 6, indicat¬ 
ing that this is a region of high stellar density compared 
to the majority of the system (see Section 3.1 for further 
discussion). 


It is clear that the emission line spectra of J2240 are 
well characterised by the superposition of multiple Gaus¬ 
sian components. The fitting residuals under the [O III] lines 
([Nil] + Ha complex) are less than I (7) per cent of the 
peak [O III] (Ha) flux. The difference in the size of the 
residuals between spectral regions is not indicative of poorer 
quality fits to the [Nil] + Ha complex, but reflects the in¬ 
crease in continuum level fluctuation with wavelength. The 
residuals under the [NH] + Ha complex are less than 2.5 
times the standard deviation of the underlying continuum. 


E 
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Ay Extinction 

Figure 3. Map of the V-band extinction {Ay) along different 
lines of sight to J2240, derived using the Balmer decrement and 
assuming an unreddened Hq;/H/ 3 ratio of 2.86. Contours trace 
the total [O HI] emission (summed across all kinematic compo¬ 
nents). The majority of the [O HI]-emitting region of the galaxy 
has Ay < 1.25 mag, with the exception of the nucleus which is 
more heavily attenuated {Ay ~ 1.75 mag). 


whereas this number can be as large as 9 under the [O III] 
lines due to their large equivalent widths. 

The quality of the fits to the [S 11] doublet indicates 
that the skylines do not compromise the accuracy of our 
emission line fitting in the red channel. However, the derived 
errors on the [S II] line fluxes are unrealistically high due to 
the presence of a skyline in the middle of the doublet. We 
therefore accept or reject the fits to the [S 11] line based on 
the signal-to-noise of the [Nil] A6583 line. 


2.5 Extinction correction 


We calculate the V-band extinction {Ay ^ in magnitudes) 
along the line of sight to each spaxel using the Ha/H/3 flux 
ratio. The hydrogen Balmer series originates from recom¬ 
bination to the n = 2 level of hydrogen. The unreddened 
intensity ratios between different lines in the series are deter¬ 
mined by their transition probabilities which vary depending 
on the temperature and density of the line-emitting gas as 
well as the shape of the ionizing radiation field. The extinc¬ 
tion along any line of sight can be calculated by comparing 
measured Balmer line ratios with their intrinsic values. 

We convert the Ha/H^d ratio in each spaxel to an Ay 
value using the [Eischera Dopit^ (|2QQ5 |) extinction curve 


with Ry = 4.5 (see Vogt, Dopita & Kewley 2013 for a 
summary of the calculation method). We adopt an unred¬ 
dened Ha/H/3 ratio of 2.86, appropriate for Case B re¬ 
combination in a nebula with an electron temperature of 
10,000 - 20,000 K and an electron density of ~I00 cm“^ 
( Osterbrock fc Eerland|2006 ). The electron temperature and 
density ranges are selected based on the findings of |Schirmer| 
et al.|(|20I3|. 


The resulting Ay map is shown in Eigurej^ Contours 
in this and all subsequent figures trace the total [O III] 
emission (summed across all kinematic components) un¬ 
less otherwise stated. The majority of the [O HI]-emitting 
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Blue Green Red 



Figure 4. Maps of the integrated continuum emission in the (left to right) blue, green and red channels. Contours trace the total [O III] 
emission. The continuum reveals the underlying stellar body of the galaxy, which is mostly rounded in shape but somewhat elongated 
towards the north eastern (right) side of the system. There is very little continuum emission in the vicinity of the bright line-emitting 
cloud, suggesting that this ionized gas cloud may be spatially distinct from the main body of the galaxy. 



Figure 5. Maps of [Nelll] A3869, Hell A4686 and [NeV] A3426 flux. Contours trace the total [OIII] emission. Black plus signs mark 
the location of the continuum peak. The [Nelll] flux has very similar spatial distribution to [OIII], which is expected given that the 
difference in their required ionization energies is only 5.84 eV. There are two peaks in the Hell emission: one coincident with the [Nelll] 
peak and one coincident with the continuum peak. The [NeV] emission, requiring the highest ionization energy, is spatially conflned to 
a small region which allows us to pinpoint the location of the central AGN. Other regions with prominent Hell and/or [NelH] emission 
are likely to be ionized by shocks. 


region of the galaxy has 0.75 < Ay < 1.25 mag, with the 
exception of the nucleus which is more heavily attenuated 
{Ay ~ 1.75 mag). The disc region is more heavily attenuated 
than the cloud, consistent with the findings of|Schirmer et aL| 
( 2 ^. 

We also calculate extinction values using an increased 
unreddened Balmer decrement of 3.1 which accounts for the 
collisional excitation of Ka in the presence of the AGN ioniz¬ 
ing radiation field ( Gaskell Ferland||1984 ). This systemat¬ 
ically decreases the calculated extinction values by 0.3 mag. 
Throughout this paper we correct line fluxes for extinction 
using both sets of extinction values and discuss the impact 
on our derived quantities. 


3 IONIZATION STRUCTURE AND SOURCES: 
BUILDING UP THE PICTURE OF J2240 

Our emission line fitting indicates that several kinematic 
components contribute significantly to the emission signa¬ 
ture of J2240. However, determining the geometry and ion¬ 
ization source of each kinematic component can be a non¬ 


trivial task, especially for systems in which the line profiles 
and dominant ionizing source(s) vary significantly across the 
galaxy. In the following section, we use integrated contin¬ 
uum maps to trace the stellar body of the galaxy and map 
the flux in high ionization energy emission lines to pinpoint 
the location of the central AGN. We then compare maps of 
the continuum and [O III] A5007 fluxes over several velocity 
slices covering the full dynamical range of the system. This 
allows us to dissect the spatial extent and relative contri¬ 
bution of individual ionizing sources as a function of both 
wavelength and position in the galaxy. 


3.1 Continuum properties 

GFHT imaging and VLT/XSHOOTER spectra of J2240 
suggest that there is a large amount of ionized gas on the 
south western side of the system which is spatially distinct 
from the main body of the galaxy (Schirmer et al. 2013). 
We investigate the location and spatial extent of the stel¬ 
lar component of J2240 by mapping the integrated contin¬ 
uum emission in each of the spectral channels (see Figure]^. 
The clear continuum peak reveals the location of the galaxy 
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centre, and is coincident with one of the two [O III] peaks 
evident from the contours. The continuum flux decreases 
smoothly with distance from the centre of the galaxy. An 
elongation is visible in the continuum emission towards the 
north eastern side of the system. A similar feature is seen in 
the [O III] emission, suggesting that the system may have 
been morphologically disturbed during recent galaxy-scale 
interactions. We note that there is very little continuum 
emission on the south western side of the system, despite 
the strength of the nebular emission lines observed in this 
region (see e.g. Figure [^. The limited extent of the stellar 
continuum emission is consistent with the hypothesis that 
the ionized gas cloud extends beyond the stellar body of the 
galaxy. 


3.2 Uncovering the AGN using Hell, [Nelll] and 
[Ne V] 

Emission line ratios extracted from the narrow line region 
of J2240 indicate that the central engine is likely to be pow¬ 
ered by an AGN (see e.g. appendix A of [Schirmer et al. 
2013). However, a single ID VLT/XSHOOTER spectrum 


across the major axis of the galaxy is insufficient to de¬ 
duce the location of the central AGN, which could be off¬ 
set from the slit axis. Our lEU data eliminate this issue 
by allowing us to investigate how the hardness of the ion¬ 
izing radiation field varies across the entire optical extent 
of the galaxy. Doubly ionized oxygen (excited by an ioniza¬ 
tion energy of 35.12 eV) is detected across a very spatially 
extended region of J2240, but probing species excited by 
larger ionization energies allows us to pinpoint the hardest 
ionizing sources in the galaxy. In particular, we investigate 
the distribution of flux in [Nelll] A3869, Hell A4686 and 
[NeV] A3426, which are excited by ionization energies of 
40.96 eV, 54.4 eV and 97.12 eV respectively. Neon is a rel¬ 
atively abundant element, produced by carbon burning in 
the late stages of the evolution of massive stars ( [Levesque Sz 
Richardson]2014| ) . The strength of the [Nelll] emission line 
correlates with the [O III] flux, and therefore [Nelll] can be 
used in place of [O HI] in AGN diagnostics ([Perez-Montero 
^rar]|2007t . However, the relatively low ionization energy 
of [Ne H] makes [Nelll] sensitive to slow and fast shocks. 
[NeV] requires a very large ionization energy and is there¬ 
fore only detectable in the immediate vicinity of an AGN 
or in material with shock velocity greater than 500 km s“^ 
( Allen et al.|2008' ). Helium is also a very abundant element, 
produced from hydrogen burning during the main sequence 
lifetime of stars. He 11 can be excited by AGN activity, fast 
shocks and young, hot stars (such as O-stars or Wolf-Rayet 
stars). Given the weakness of the stellar continuum emis¬ 
sion in J2240, it is unlikely that active star-format ion con¬ 
tributes significantly to the Hell emission. Hell is excited 
by an ionization energy between that of [NelH] and [NeV], 
and therefore we expect the Hell flux distribution to be 
intermediate between that of the two neon lines. 

We show maps of the total [Nelll], Hell and [NeV] 
emission in Eigurej^ (It is not necessary to break the flux 
into components because these emission lines only originate 
from the most energetic ionizing sources.) The spatial distri¬ 
bution of the [NelH] emission is very similar to that of the 
total [O HI] emission within the main body of the galaxy, 
suggesting that the ionization parameter is high enough 


across the galaxy to excite both lines equally. However, the 
[NelH] emission in the cloud is concentrated around the 
northern and southern [O HI] peaks, suggesting that these 
regions may be shock-excited. 

The Hell flux is more concentrated towards the centre 
of the galaxy than the [Ne HI] flux, and has a double-peaked 
distribution. One of the peaks is coincident with the [NelH] 
peak, whilst the second peak is coincident with the contin¬ 
uum peak. We also see evidence for Hell emission in the 
southern [O HI] peak of the cloud, strengthening the case 
for shock excitation in this region. The [NeV] emission is 
very spatially concentrated around the continuum peak. The 
presence of a very hard ionizing radiation field at the centre 
of the galaxy is consistent with AGN activity. The second 
Hell peak is invisible in [NeV], and is therefore likely to be 
powered by fast shocks. There is no evidence for a second 
active galactic nucleus from this data. 


3.3 [O III] A5007 velocity channel maps 

We probe the spatial and spectral variation in the emis¬ 
sion line morphology of J2240 by constructing flux maps 
for nine velocity channels across the [O III] A5007 emission 
line profile (see Eigure[^. Maps of the continuum emission 
across the same channels have been included below each 
emission line map for comparison. No signal-to-noise ratio 
cuts have been applied. Gontours on the emission line maps 
trace the total [O HI] emission, whilst contours on the con¬ 
tinuum maps trace the integrated continuum emission from 
the red channel. The colorbar is scaled independently in each 
map and velocity slice to span the full range of flux levels, 
from zero flux in black to maximum flux in red. 

It is clear that the morphology of the continuum emis¬ 
sion does not evolve significantly as a function of wavelength. 
This reflects the fact that the continuum emission can be ex¬ 
plained by a single, kinematically simple emission source - 
the underlying stellar population of the galaxy. In contrast, 
the morphology of the [O HI] emission evolves dramatically 
between velocity slices. The highest velocity gas (slices 1, 2, 
8 and 9) is primarily aligned along the major axis of the con¬ 
tinuum emission, suggesting that it may be associated with 
the underlying stellar population of the galaxy (i.e. tracing 
the gas disc of the galaxy). However, the central velocity 
slices are dominated by line emission from the gas cloud 
which is misaligned by ~ 45° compared to the continuum 
axis. The emission from the cloud is clumpy and morphologi¬ 
cally irregular, suggesting that there may be significant vari¬ 
ation in the ionization parameter and/or gas density within 
the cloud. 

The evolution in emission morphology as a function of 
wavelength suggests that a) the rotation curve of the galaxy 
disc has a larger velocity amplitude than the rotation curve 
of the cloud structure, and b) the cloud component has a 
much larger [O III] luminosity than the south western-most 
region of the disc, causing the former to be more promi¬ 
nent in the central velocity slices. Slices four and five have 
the largest [O HI] fluxes, and both are dominated by [O III] 
emission from the continuum peak, consistent with our hy¬ 
pothesis that this is the location of the central AGN. 

The channel maps reveal several kinematically broad 
regions which emit strongly across several velocity channels. 
The elongation on the north eastern side of the galaxy is 
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Figure 6. Flux maps for nine velocity channels across the [OIII] A5007 emission profile, with continuum maps extracted from the same 
channels for comparison. Contours on the emission line maps trace the total [O III] emission, and contours on the continuum maps 
trace the integrated continuum emission in the red channel. It is clear that the morphology of the continuum emission does not evolve 
significantly as a function of wavelength. In contrast, the morphology of the [O III] emission evolves dramatically between velocity slices. 
The highest velocity gas (slices I, 2, 8 and 9) is primarily aligned along the major axis of the continuum emission suggesting that it may 
be associated with the galaxy’s underlying stellar population. However, the central velocity slices are dominated by line emission from 
the gas cloud which is misaligned by ~ 45° compared to the continuum axis. 
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Figure 7. Spectra centred around the [O III] A5007 and [O III] A4959 emission lines, highlighting the superposition of two narrow 
components along various edges of the galaxy disc. The second component is associated with the galaxy cloud. The velocity offset is due 
to the differing velocity amplitudes of the cloud and disc rotation curves. 


seen in the last four velocity slices, the emission line peak 
between the disc and the cloud is seen prominently in the 
first three velocity slices, and the northern [O III] peak in the 
cloud is observed in the last five velocity slices. This suggests 
that there may be several localised regions of energetic, high 
temperature gas in J2240. 

Our analysis suggests that at least three ionizing sources 
contribute significantly to the line emission of J2240 - the 
gas disc, which extends along the major axis of the galaxy’s 
continuum emission, the gas cloud, which is ionized along 
a distinctly misaligned axis, and more energetic ionizing 
source(s) which increase the turbulent energy of the gas. 
This is consistent with the line profiles and emission line fits 
presented in Figure We expect the two narrow compo¬ 
nents to be superimposed at the edge of the gas disc, and 
observe double-peaked line profiles in several locations (see 
Figure]^. The two narrow components are offset in veloc¬ 
ity due to the different velocity amplitudes of the disc and 
cloud rotation curves, but appear to have very similar veloc¬ 
ity dispersions. Double peaked line profiles are not visible in 
the centre of the galaxy because emission from the galaxy 
component is much stronger than emission from cloud com¬ 
ponent in this region. 


our analysis of the [O III] channel maps. Three kinematic 
components are required in the vicinity of the elongation, at 
the [Nelll] peak, in the region between the cloud and the 
main body of the galaxy and at the northern [O III] peak 
of the cloud. Our analysis of the channel maps has revealed 
that 3/4 of these regions have large velocity dispersions, and 
the [Nelll] peak is detected in [Nelll] and Hell, suggesting 
that energy injection from shocks is likely to drive the gas 
turbulence in all four regions. 

The subsequent sections of this paper will delve into the 
physical properties of each of the individual kinematic com¬ 
ponents, making it imperative that the components are num¬ 
bered consistently across the galaxy according to their phys¬ 
ical origin. LZIFU automatically numbers the components in 
order of increasing velocity dispersion within each individual 
spaxel. This guarantees that the broad components identi¬ 
fied in the channel maps will always appear in component 
3. The two narrow components have similar velocity disper¬ 
sions, but the velocity dispersion in the disc is consistently 
larger than that in the cloud, ensuring that component 1 
will be dominated by emission from the ionized cloud and 
component 2 will trace the rotation curve of the main body 
of the galaxy as well as some of the turbulent gas in the 
cloud. 


3.4 Component assignment 

Figure shows the number of statistically significant kine¬ 
matic components in each spaxel of J2240 (determined as de¬ 
scribed in Section 2.4). Spaxels requiring one, two and three 
components are coloured light grey, dark grey and black re¬ 
spectively. Spaxels for which the emission line fitting did not 
converge are not coloured. At least two kinematic compo¬ 
nents are required to reproduce the observed spectra across 
the majority of the [O III]-emitting region, as expected from 


4 3D SPECTROSCOPIC ANALYSIS 
4.1 [O III] morphology 

Figure shows [O III] flux maps for each of the individual 
kinematic components. The flux range spanned by the col- 
orbar is fixed across all three maps to indicate the relative 
contribution of different kinematic components and physical 
regions to the overall [O III] emission of the system. Compo- 















Figure 8. Map showing the number of statistically significant 
kinematic components in each spaxel of J2240, as determined by 
comparing the reduced values of the 1, 2 and 3 component fits. 
Spaxels requiring 1, 2 and 3 components are coloured light grey, 
dark grey and black respectively. Spaxels for which the emission 
line fitting did not converge are not coloured. Contours trace the 
total [O III] emission. At least two kinematic components are re¬ 
quired to reproduce the observed spectra across the majority of 
the [O Illj-emitting region, and several regions which are likely to 
be shock-excited require three components. 


nent 1 clearly traces emission from the ionized gas cloud, as 
expected. There is a large amount of diffuse gas with rela¬ 
tively weak [O III] emission, and the regions of most intense 
[O III] flux appear to be strongly associated with the mis¬ 
aligned emission axis identified from the channel maps. This 
emission axis intersects with the main body of the galaxy 
at the location of the AGN, suggesting that we may be ob¬ 
serving an AGN ionization cone (indicated by the red lines 
in Figure]^. 

The majority of the [GUI] emission in component 2 
originates from the main body of the galaxy. The spatial 
distribution of [O III] flux in this component is very similar 
to the total [O III] flux distribution traced by the contours, 
indicating that the second component is the dominant source 
of [O III] flux across the majority of J2240. The [O III] emis¬ 
sion peaks strongly at the location of the central AGN. We 
also clearly recover the northern and southern [GUI] peaks 
in the cloud, which have both been identified as possible 
shocked regions. 

The [G III] emission in component 3 originates from 
a similar region of the system to component 2, but with 
significantly different morphology. The [Nelll] peak is the 
dominant source of [GUI] emission in component 3, and 
has a similar peak flux to the AGN emission in component 2. 
The strong [G III] emission from the [Ne III] peak appears to 
drive the strong warp in the [GUI] contours. The elongation 
on the north eastern side of the galaxy is clearly revealed in 
this component, and appears to be physically distinct from 
the remainder of the galaxy. We also observe a region of 
weaker, diffuse emission between the north western edge of 
the galaxy disc and the northern [GUI] peak of the cloud. 
This feature is not observed in either the channel maps or 
the high ionization line maps, and is only revealed when the 
much stronger narrow components are removed. 
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4.2 Velocity field and velocity dispersion 


Figure 10 shows the velocity and velocity dispersion (a) 
fields of each of the kinematic components. We also plot 
position-velocity and position-cr diagrams extracted along 
the major axis of the galaxy down the right hand column 
of the figure. Positions are calculated with respect to the 
central AGN and components 1, 2 and 3 shown in red, blue 
and green respectively. 

The velocity field of component 2 reveals a clear rota¬ 
tion curve over the disc region of the galaxy with an am¬ 
plitude of ~300 kms“^. The velocity field of component 1 
is also dominated by rotation, but the posit ion-velocity dia¬ 
gram reveals that its rotational velocity is smaller than that 
of component 2, as expected from our analysis of the chan¬ 
nel maps. The co-rotation of component 1 and component 2 
suggests that the cloud material which is superimposed over 
the main body is enveloping it and is therefore morpholog¬ 
ically linked to the main body of the galaxy (as discussed 
by Schirmer et al. 2013). The remainder of the cloud ma¬ 


terial generally lies at the systemic velocity of the galaxy, 
although some shocked regions in component 2 have more 
turbulent emission components with anomalous velocities. 
The emission from the northern [G III] peak in the cloud 
is redshifted by about 200 kms“^, and there is some evi¬ 
dence for blueshifted material in the north eastern region of 
the cloud. The average velocity dispersions of components 
1 and 2 are 100 kms“^ and 150 kms“^ respectively, indi¬ 
cating that the gas is globally disturbed despite the ordered 
rotation we observe. 

Gomponent 3 traces the most turbulent gas in the sys¬ 
tem. The average velocity dispersion is 300-400 kms“^, but 
two clumps of gas (in the vicinity of the elongation and 
between the two [G III] peaks in the cloud) have anoma¬ 
lously high velocity dispersions of 600-800 kms“^. These 
distinct clumps of high velocity dispersion gas are evidence 
for localised energetic phenomena such as shocks or out¬ 
flows, whilst the majority of the underlying broad compo¬ 
nent is likely to be energized by the central AGN. There 
is some evidence for rotation in component 3; however, the 
kinematics to the north east of the central AGN are very 
irregular. A strongly blue-shifted clump of gas is visible just 
to the south of the [Nelll] peak, and gas in the vicinity of 
the elongation appears strongly redshifted. 


4.3 Diagnostic diagrams 

The [NII]/Ho and [SII]/Ho vs. [GIII]/H/3 diagnostic dia¬ 
grams are commonly used to separate star formation from 
shock excita tion and AGN activity ([Baldwin, Phillips Ter- 


levich 1981 Veilleux & Gsterbrock 1987| ). The diagnostic 


line ratios are virtually insensitive to reddening due to the 
small wavelength separation of the line pairs ([NII] A6583 & 
Ra A6563, [S II] AA 6717,6731 k Ra A6563, [GUI] A5007 k 
H/3 A4861). The [GIII]/H/3 ratio is sensitive to the electron 
temperature and ionization parameter of the line-emitting 
gas, and the [Nil]/Ho ratio is a tracer of metallicity due 
to the secondary nucleosynthetic production of nitrogen at 
high oxygen abundances (see e.g. [Kewley et al.|2013| ). 

Figure shows [Nil]/Ho (left-hand column) and 
[SH]/Ho (right hand column) vs. [GIH]/H/3 diagnostic di¬ 
agrams with line ratios extracted from the spectrum of each 
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Figure 9. [O III] flux maps for each of the three line-emitting components. Contours trace the total [O III] emission. The kinematic 
components have different [O III] morphologies, suggesting that they have different physical origins. Component I traces emission from 
the ionized gas cloud, component 2 traces the rotation curve of the main body of the galaxy and component 3 traces turbulent gas 
energized by fast shocks. Red lines overplotted on the component I map trace the edges of the possible ACN ionization cone. 
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Figure 10. Maps of the velocity held and velocity dispersion (cr) for each of the three line-emitting components, and position-velocity 
and position-cr diagrams extracted along the major axis of the galaxy. Contours trace the total [O III] emission. In the two right hand 
panels, positions are with respect to the central AGN and components I, 2 and 3 are shown in red, blue and green respectively. All three 
components show clear rotation signatures, but the gas is globally disturbed. The velocity dispersion reaches 600-800 kms“^ in some 
regions of component 3, indicating energy injection from fast shocks and outffowing material. 


kinematic component in each spaxel of J2240. Components 
1, 2 and 3 are shown in red, blue and green respectively. The 
Kewley et al. |( |20oT] i theoretical upper bound to pure star- 
formation (solid lines) separates star-formation (bottom left 
of diagrams with lower [NII]/Ha, [SII]/Ha and [OIII]/H/3 
ratios) from harder ionizing sources (top right of diagrams 
with higher [NII]/Ha, [SII]/Ha and [OIII]/H/3 ratios). All 


of the observed line ratios for each of the kinematic com¬ 
ponents lie strongly in the AGN region of the diagnostic 
diagrams. 

We observe that the [O III]/Hy5 ratios are inversely cor¬ 
related with both the [NII]/Ha and [SII]/Hq; ratios. We 
explore possible origins of this ionization signature by plot¬ 
ting Dusty AGN (top panel; [Groves, Dopita Sutherland 


2QQ4| , fast shock (mi ddle panel; [Allen et al.|2Q08| ) and slow 


shock (bottom panel; [Rich et al.|2QlQ[ ) models over the data, 
shown in brown. The selected AGN models are calculated 
for a metallicity of 12 + log(0/H) = 8.99 (within the 0.9- 


2Zq range derived by Schirmer et al.|2013 ) and an electron 
density of rie — 100 cm“^ (between the <50 cm“^ density 
of the cloud and the 150-650 cm“^ density of the disc in 
J2240; [Schirmer et al.[[2013| . Solid lines connect points of 
constant power-law index (increasing from T = 1.2 - 2 as the 
[Nil]/Ha ratio decreases), and the ionization parameter in¬ 
creases from log U = -4 to log = 0 as the [O III]/H/3 ratios 
increase along each line. It is clear that ionization parameter 
variations are able to produce the observed anti-correlation 
in the line ratios (see Section [P^ for further discussion of the 
ionization parameter). Variations in the metallicity may ac¬ 
count for the spread in [NII]/Ha ratios at a given [O III]/H/3 
ratio. 

The fast shock models are calculated for a metallicity 
of 12 + log(0/H) = 8.69 and a pre-shock hydrogen density 
of uh = 100 cm“^ (consistent with the metallicity and den¬ 
sity ranges derived by Schirmer et al.|2013 ). Solid lines con¬ 
nect points of constant magnetic field strength (ranging from 
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Figure 11. [NII]/H a and [SII]/Hq; vs. [0III]/H^ diagnostic diagrams with line ratios extracted for each component of each spaxel. 
Red, blue and green points are extracted from c omponents 1, 2 and 3 respectively, a nd brown points have been extracted from model 
grids. The top panel shows Dusty AGN models ( [Groves, Dopita &: Sutherland||2004| with 12 -h log(0/H) = 8.99 and Ue = 100 cm“^. 
Solid lines connect points of constant power law index, ranging from F = 1.2-2. The middle panel shows fast shock models ( [Allen et al.j 
[2008[ > with 12 -h log(0/H) = 8.69 and precursor density of nn = 100 cm“^. Solid lines connect points of constant magnetic field strength, 
going from lO^'^-lOO p,G. The bottom panel shows slow shock models ( [Rich et al.[[2010| >. Solid lines connect points of consant shock 
velocity, ranging from 100-200 kms“^. The Dusty AGN models indicate that the anti-correlation between the [Nil]/Ha and [OIII]/H^ 
ratios can be explained by ionization parameter variaions in the line-emitting gas. The fast shock models also coincide with a significant 
fraction of the data points, indicating that fast shocks could be the principal ionizing source in some regions of the galaxy. 


10“^-100 /xG). The [OIII]/H/3 ratios increase as the shock 
velocity increases from 300-650 kms“^ and then decrease 
again slightly as the shock velocity approaches 1000 kms“^. 
The fast shock models lie in the same region of the diagnos¬ 
tic diagram as data points with log([0 III]/Hy5) < 1. This 
suggests that fast shocks could contribute significantly to 
the ionization signature in some regions of J2240, consistent 
with our identification of several prominent shocked regions 
throughout the system. 


The slow shock models are calculated for a range of 
metallicities (8.39 < 12 + log(0/H) < 9.29) and shock ve¬ 
locities (100-200 kms“^). Solid lines connect points of con¬ 
stant shock velocity, increasing towards larger [Nil]/Ha ra¬ 
tios, and the metallicity increases towards larger [Nil]/Ha 
ratios along each line. It is clear that slow shocks are not 
powerful enough to produce the large [OHI]/H/3 ratios we 
observe, indicating that even if they are present in J2240 
they are not the dominant ionizing source. 


4.4 Maps of diagnostic line ratios 


Mapping diagnostic line ratios across galaxy systems reveals 
spatial variations in the ionization state of the gas. Figure 


12 shows maps of the [NII]/Ha, [SII]/Ha and [OHI]/H/3 


line ratios for each of the kinematic components. The av¬ 
erage [Nil]/Ha ratios increase from component 1 to com¬ 
ponent 3 (-0.236, -0.211, -0.120), and the same behaviour is 
seen in the [SII]/Ha ratios (-0.271, -0.234, -0.217). The av¬ 
erage [OHI]/H/3 ratios decrease from component 1 to com¬ 
ponent 3 (1.003, 0.955, 0.882), demonstrating the observed 
anti-correlation between the [O IH]/H/3 ratios and [NII]/Ha 
and [SII]/Ha ratios. 


The component 1 [OIH]/H/3 map has a similar mor¬ 
phology to the [O HI] flux map in Figure There is very 
little variation in the [O HI]/H/3 ratio over the [O HI]-bright 
region of the galaxy (0.9 < log([0 HI]/H/3) < 1.1), suggest¬ 
ing that the temperature and ionization parameter of the 
line-emitting gas are relatively uniform. The [Nil]/Ha and 
[SH]/Ha ratios peak at the elongation and various regions 
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Figure 12. Maps of the [NIIJ/Ha, [SII]/Hq; and [OIIIJ/H/3 ratios for each of the individual components. Contours trace the total 
[O III] emission. Black plus signs indicate the location of the central AGN. The average [Nil]/Ha and [SII]/Ha ratios increase from 
component I to component 3, whilst the average [OIII]/H/3 ratio decreases. The largest [Nil]/Ha and [SII]/Ha ratios are found on the 
north eastern side of the system, consistent with our hypothesis that fast shocks are prominent in this region. 


in the cloud, suggesting that these regions may have high 
gas densities and/or significant shock excitation. 

There is clear spatial variation in the [OIII]/H/3 ratios 
of component 2. The largest [OIII]/H/3 ratios are observed 
in the vicinity of the central AGN and along the ionization 
cone. The [OIII]/H/3 ratios decrease steadily with distance 
from the axis of the ionization cone, driven by a correspond¬ 
ing decrease in the number of ionizing photons. In contrast, 
the [Nil]/Ha and [SII]/Ha ratios increase with distance 
from the axis of the ionization cone, consistent with our 
hypothesis that the observed anti-correlation between the 
diagnostic line ratios is driven by ionization parameter vari¬ 
ations. 

We do not observe a strong relationship between the 
component 3 line ratios and position relative to the AGN 
ionization cone. The largest [Nil]/Ha and [SII]/Ha ratios 
are found towards the north eastern side of the galaxy, in 
the vicinity of the elongation and the [Nelll] peak. These 
large line ratios support our hypothesis that the gas in this 
region of the galaxy is primarily powered by fast shocks. 


4.5 Ionization parameter 

The ionization parameter q is defined as the ratio between 
the flux of ionizing photons and the local hydrogen number 
density. This ratio is a measure of the maximum velocity at 
which the radiation field can drive an ionization front. The 
more commonly used dimensionless ionization parameter U 
is related to q by U = q/c (where c is the speed of light). 


The combination of ionization parameter and electron den¬ 
sity can be used to derive the photon flux (luminosity) re¬ 
quired to ionize a particular region. In the following sections 
we explore the ionization parameter, electron temperature 
and electron density of the line emitting gas in J2240. De¬ 
tailed calculations of the ionizing luminosity of the EELR 
are deferred to a future paper. 


The ionization parameter in HII regions can be calcu¬ 
lated from the [OIII]/[OII] ratio given an estimate of the 
abundance of the line emitting gas ( Kewley Dopita|2QQ2 ). 
However, collisional excitation in the presence of an AGN in¬ 
creases the [OIII]/[OII] ratio for a given ionization param¬ 
eter, making the conversion between [OIII]/[OII] ratio and 
ionization parameter difficult. Mapping the [O III]/[0 II] ra¬ 
tio can provide some insight into the spatial variation in the 
ionization parameter, although calculation of the absolute 
ionization parameter values is very uncertain. 


We show maps of the [OIII]/[OII] ratios for each of 
the kinematic components in Eigure The [OIII]/[OII] 
ratio in component 1 peaks at the location of the [Nelll] 
peak. To the south west of the AGN, the [OIII]/[OII] ra¬ 
tio is largest along the axis of the AGN ionization cone and 
decreases steadily with distance off the axis. This indicates 
that the number of ionizing photons decreases with increas¬ 
ing opening angle, consistent with expectations for an ion¬ 
ization cone. 


Gomponent 2 shows very similar [O III]/[0 II] ratios to 
component 1, consistent with our hypothesis that both com¬ 
ponents are ionized primarily by the central AGN. The peak 
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Figure 13. Maps of the [OIII]/[OII] ratios for each of the kinematic components. Contours trace the total [O III] emission. The 
[O III]/[Oil] ratio can be used to calculate the ionization parameter in HII regions but is enhanced by collisional excitation in the 
presence of an AGN. In components 1 and 2 the ionization parameter decreases with distance from the axis of the AGN ionization cone, 
as expected. The average [O III]/[Oil] ratios in component 3 are significantly lower than those of the first two components, reflecting 
the prominence of shock excitation in ionizing the most turbulent gas in the system. 


[O III]/[Oil] ratios in component 2 occur in the vicinity of 
the central AGN, but large errors on the [O II] flux values 
prevent robust calculation of the ratio over a significant frac¬ 
tion of the nuclear region. The average [OIII]/[OII] ratios 
in component 3 are significantly lower than those of the first 
two components. This is consistent with a larger contribu¬ 
tion from shocks, which can have lower [OIII]/[OII] ratios 
than Seyfert 2 nuclei (e.g. Dopita Sutherland||1995 ). 

The Dusty AGN models shown in Figure mr can be 
used as a quantitative probe of the ionization parameter val¬ 
ues across J2240. The line ratios of AGN-dominated spaxels 
(log([N II]/H(a) < 0) are best reproduced by models with 
\ogU < —1.3 at a = 1.7 or —2.6 < \ogU < —1.3 at a = 1.4. 
The true a and log U values are likely to lie between these 
two estimates. The ionization parameter values we derive 
are larger on average than the values calculated by |Schirmer| 


et al. (2013). This is likely because we assume a metallicity 


of ~ 2 ^ 0 , where Schirmer et al. (2013) base their calcula¬ 


tions on the Kewley & Dopita (2002) ionization parameter 
calibration assuming a metallicity of ~ 0.5 Zq. 


4.6 Electron temperature 

The temperature and density of electrons are also power¬ 
ful probes of the principal ionization source(s) of the line- 
emitting gas. Gompression and heating of gas along shock 
fronts can produce electron temperatures above 20,000 K 
(e.g. Evans et al. 1999), compared to an average tempera¬ 
ture of ~10,000-15,000 K in AGN narrow line regions (e.g. 
[Oster brock fc Ferland|20(36 ). This provides a means to sep¬ 
arate ionization due to fast shocks from ionization due to 
AGN activity. 

The electron temperature determines the relative pop¬ 
ulation rate for different excited states of the same atomic 
species. The [GUI] A4363 / ([GUI] A5007 + [GUI] A4959) 
ratio provides an estimate of the total cooling due to 
oxygen in the medium ionization zone of nebulae, and is 
therefore commonly used as an estimate of the electron 
temperature (Gsterbrock & Ferland 20Q6| ). [GUI] A5007 


and [GUI] A4959 have excitation temperatures of 29,130 K 
whereas [G III] A4363 is an auroral line with an excitation 
temperature of 62,094 K. Therefore, the intensity ratio of 
the auroral line to the two oxygen strong lines increases 
with temperature. 

The relationship between the measured [GUI] line ra¬ 
tio and the electron kinetic temperature depends on both 
the electron density and the intrinsic electron tempera¬ 
ture distribution. The rate of collisional excitation increases 
with density, and therefore the [GUI] ratio will be larger 
in high density nebulae than in low density nebulae at 
a given electron temperature. However, at densities below 
He = 10 ^cm“^ this effect is negligible and the [G III] ra¬ 
tio becomes independent of density ( Nicholls et al.|[2013| ). 
Schirmer et al. (2013) find that the electron density in J2240 
rarely exceeds 650cm“^, indicating that electron density 
does not need to be considered in our electron temperature 
calculations. 

Traditionally, electron temperature calibrations have 
been based on the assumption that electron energies follow 
a Maxwell-Boltzmann distribution. However, | Nicholls, Do-| 
pita & Sutherland ( 2012 ) and Nicholls et al. ( |2013| ) showed 
that a At distribution is more appropriate for H H regions in 
the local universe. The k value controls the strength of the 
high energy electron tail. A /t value of infinity corresponds to 
the Maxwell-Boltzmann distribution, and the fraction of hot 
electrons increases as the k value decreases. A larger fraction 
of hot electrons will increase the expected [G III] ratio rela¬ 
tive to the Maxwell-Boltzmann estimate for a given kinetic 
temperature, causing calibrations based on the Maxwell- 
Boltzmann distribution to over-estimate the true nebular 


electron temperature. Nicholls, Dopita & Sutherland ( 2012 ) 


show that correcting the electron temperature calibrations 
for the increased hot electron fraction appropriate to a k- 
distribution resolves the long standing discrepancy between 
electron temperatures calculated using direct and strong line 


ratio methods. The correction is as follows (Nicholls et al 


2013 see their Equation 36): 
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Figure 14. Maps of (top) electron temperature, (middle) error on the electron temperature and (bottom) electron density for each of the 
three kinematic components. Contours trace the total [O III] emission. Component 3 has the highest average electron temperature. The 
[Nelll] peak, continuum peak and elongation all show electron temperatures between 20,000-25,000 K, supporting our hypothesis that 
shock excitation is significant in these regions. Much of the diffuse gas in the cloud is in the low density limit, however, we are able to 
derive density measurements for the [O Illj-bright regions along the AGN ionization cone, suggesting that material is being compressed 
along the ionization cone. The density reaches several hundreds to thousands of cm“^ in the hottest, most energetic regions of the system. 


Tkin = a{Hl)Te,MB + b{B) 


( 1 ) 


where Te,MB is the electron temperature estimate ac¬ 
cording to the original Maxwell-Boltzmann calibration of 


Osterbrock & Ferland (2006) and a(/^) and 6(/^) are coeffi¬ 
cients which depend on n. We adopt a = 14192.273 and b = 
2413.5652, calculated using Equations 37 and 38 of |NichoIIs| 
et al. (2013) with b = 20. The average b values of AGN nar¬ 
row line regions are expected to lie between b = 10-20. Zones 
of very low b values (adjacent to regions of very high exci¬ 
tation) mix with equilibrium zones, producing intermediate 
average b values (David Nicholls, private communication). 
Schirmer et al. I poT^ analysed the electron tempera¬ 


tures in J2240 using the [O III] Te proxy, but due to the 
weakness of the [O III] A4363 auroral line in their data they 
were only able to derive globally averaged temperatures for 
the cloud and main body of the galaxy. They found the main 
body of the galaxy to have a higher average temperature 
than the cloud (17,100 K compared to 14,700 K). However, 
such a globally averaged analysis masks any smaller scale 


temperature variations which would trace the localized im¬ 
pact of shock fronts on the energetics of the gas. Our CMOS 
data are much deeper and allow us to not only map the elec¬ 
tron temperature across the galaxy, but also to do this on a 
component by component basis. 


Unfortunately, we are unable to detect [O III] A4363 at 
the 3cr level in every component of every spaxel. We sac¬ 
rifice a small amount of our spatial resolution by binning 
spaxels to achieve a target signal-to-noise ratio of 3. Tradi¬ 
tional binning techniques decrease the spatial resolution over 
all regions of the galaxy uniformly, regardless of variations 
in the initial signal-to-noise ratio across the galaxy. Voronoi 
binning is an adaptive spatial binning technique which pro¬ 
duces the optimal binning solution by tessellating bins of 
varying size ( Cappellari Copin|2QQ3 ). Spaxels whose orig¬ 
inal signal-to-noise ratios exceed the target are not binned, 
and the bin sizes increase towards lower signal-to-noise re¬ 
gions. 


Figure shows electron temperature maps for each 
of the kinematic components. Component 3 has the high¬ 
est average electron temperature, consistent with our iden- 
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tification of turbulent, shock- and AGN-ionized gas in this 
component. The [Nelll] peak, the elongation and the con¬ 
tinuum peak all have temperatures between 20,000 K and 
25,000 K, consistent with our hypothesis that these regions 
are excited by a combination of shock excitation and AGN 
activity. 

The first two components, by contrast, are significantly 
cooler. The gas in component 1 has an average tempera¬ 
ture of 14,700 K, and the gas in component 2 has an aver¬ 
age temperature of 12,900 K, consistent with expectations 
for AGN narrow line region gas. Grouping the temperature 
measurements by physical region, we find that the average 
electron temperatures in the cloud and the main body of the 
galaxy are 13,700 K and 13,900 K respectively. We measure 
an 800 K temperature difference between the two compo¬ 
nents, which is much smaller than the 2,400 K difference 


found by Schirmer et al. (2013). However, if we include the 


more energetic third component in our calculations then the 
temperature of the main body of the galaxy rises to 15,500 K 
and the temperature difference between the regions rises to 
1,600 K. 

The temperatures we derive for each of the compo¬ 
nents are systematically ~ 1,500 K lower than the mea¬ 
surements of Schirmer et al. (2013). This is most likely 


to arise due to a combination of our multiple-component 
parametric emission-line fitting and our assumption of k- 
distributed electrons. We check the consistency of our cal¬ 
culations by measuring the electron temperatures from the 
total [O III] A5007, [O III] A4959 and [O III] A4363 fluxes in 
each spaxel (integrated across all three kinematic compo¬ 
nents), assuming a standard Maxwell-Boltzmann electron 
temperature distribution. We find that the cloud and the 
main body of the galaxy have average electron tempera¬ 
tures of 15,600 K and 16,900 K respectively, much closer to 


the values published by Schirmer et al. (2013). 


Gorrecting for extinction based on an unreddened 
Balmer decrement of 3.1 results in systematically lower 
electron temperatures across the galaxy. The component 
1, 2 and 3 temperatures decrease by an average of 250 K 
(1.6 ± 0.4 per cent), 195 K (1.5 ± 0.3 per cent) and 480 K 
(2.2 ± 0.8 per cent) respectively. These differences are a 
very small percentage of the total electron temperatures and 
therefore the choice of unreddened Balmer decrement does 
not have a significant impact on our electron temperature 
measurements. 


fc Ferland||2006| for a more detailed explanation). The [S II] 
ratio {R[sii]) also depends on temperature, according to the 
equation 


R[sii] — a 


1 bx 
1 + cx 


( 2 ) 


where a, b and c are constants and x = 10~^. 
We use the electron temperatures derived in the previous 
section as inputs to our electron density calculations. Be¬ 
cause the [SII] doublet is not sensitive to electron den¬ 
sity for He < 100 cm“^, we present only an upper limit of 
Ue < 100 cm“^ for spaxels in which the estimated electron 
density is ^ 100 cm“^. 

We show maps of the electron density for each of the 
kinematic components in the bottom row of Figure We 
cannot calculate errors on the electron densities for any of 
the components due to the presence of skyline residuals in 
the region of the [S II] line profile (see discussion in Section 

The density maps highlight similar features to the elec¬ 
tron temperatures maps. The average electron density in 
component 1 is 220 cm“^, but this is a strict upper limit 
given that 60 per cent of the spaxels have only an upper 
limit on the electron density. The majority of the diffuse gas 
in the cloud is in this low density limit. However, we derive 
density measurements for the northern and southern [O III] 
peaks in the cloud, a clump of gas between the cloud and 
the main body of the galaxy, the [Nelll] peak and the elon¬ 
gation on the north eastern side of the system, which are 
all expected to be dense regions given their temperatures, 
velocities and line ratios. 

The average electron density in component 2 is also 
220 cm“^, but we again find that a large fraction of the 
spaxels (65 per cent) are in the low density limit and there¬ 
fore the actual average density will be lower. The density 
clearly peaks along the axis of the AGN ionization cone, 
reaching Ue ~ 700 cm“^. This is consistent with the findings 


of Schirmer et al. (2013) and suggests that AGN activity is 


contributing to the compression of material. 

The average electron density in component 3 is much 
larger, reaching into the thousands due to the presence of 
some highly compressed gas. We observe that the densest gas 
in J2240 is also the hottest and the most energetic; indicat¬ 
ing that all of the tracers and diagnostics we have explored 
paint a coherent picture of this system. 


4.7 Electron density 

The [S H] A6717/[S II]A6731 ratio is a probe of the electron 
density in the line-emitting gas. The two [S II] lines result 
from transitions with similar ionization energies, however, 
the ratio between the maximum occupancies of the transi¬ 
tions is 1.5 and the ratio between their lifetimes is ~ Iii 
low density nebulae (ue < 100 cm“^), collisional excitation 
is very rare and the ratio between the transition strengths 
is approximately equivalent to the ratio of their maximum 
occupancies. At very high densities (ne > 10, 000 cm“^), the 
ratio of the transition strengths is approximately equal to 
the ratio of the maximum occupancies divided by the ratio of 
the lifetimes (0.44). However, at intermediate densities the 
probability of the [S H]A6717 transition increases relative to 
the probability of the [SII]A6731 transition (see |Osterbrock 


5 DISCUSSION 

5.1 Piecing together the puzzle: the 3D ionization 
structure of J2240 

Our analysis of the 3D kinematic and ionization structure 
of J2240 has confirmed that this system is very complex. 
Two narrow components and an underlying broad compo¬ 
nent are ubiquitous across the majority of the system, in¬ 
dicating that at least two physical gas components and at 
least two distinct ionizing sources contribute significantly to 
the emission signature of J2240. In the following section, we 
synthesise our results from Sections and ^to piece together 
a picture of J2240. 

Our discussion is complemented by an interactive 3D 
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Figure 15. Interactive 3D contour plot, produced in .wrl format using the Python package Mayavi2, converted to .u3d format using 
PDF3D ReportGen and embedded in the PDF file using the DTgX movie15 package. By default, a 2D place-holder figure is displayed, 
but interactivity can be enabled in supported PDF readers by clicking on the figure. The interactivity allows the reader to rotate, pan 
and zoom in/out of the 3D plot. The plot shows iso-intensity surfaces for an 800 kms“^ slice of the reduced data cube extracted around 
the [O III] A5007 emission line. The red axis is velocity, going from -420 kms“^ to 380 kms“^, the green axis is Y position and the blue 
axis is X position. There are three pre-set views: position-position, velocity - Y position and X position - velocity, which can be selected 
from the “Views” drop-down menu. Two spheres embedded in the isosurfaces mark the position of the central AGN and the hotspot 
energized by interaction of the radio jet with the interstellar medium. These spheres are labelled but the text is not visible from all 
viewing angles. 



contour plot designed to assist the reader in their un¬ 
derstanding of the velocity structure of J2240 (see Fig¬ 
ure [15 ). The plot was prod uced in .wrl format using the 


Python module Mayavi^ ( jRamachandran Varoquaux 


2011), converted to .u3d format using the proprietary soft¬ 
ware PDF3D ReportGeiJ^ and embedded in the PDF file 


Cfpa 


using the DT^]X MOVIE 1^ package. By default, a 2D place¬ 
holder figure is displayed, but interactivity can be enabled in 
supported PDF readers by clicking on the figure. The inter¬ 
activity allows the reader to rotate, pan and zoom in/out of 
the 3D plot. The 3D plot shows iso-intensity surfaces for an 
800 kms“^ slice of the reduced data cube extracted around 
the [O III] A5007 emission line. The red axis is velocity, going 
from -420 kms“^ to 380 kms“^, the green axis is Y position 
and the blue axis is X position. There are three pre-set views: 
position-position, velocity - Y position and X position - ve¬ 
locity, which can be selected from the “Views” drop-down 
menu. Two spheres embedded in the isosurfaces mark the 
position of the central AGN and the hotspot (discussed at 
the end of this section). These spheres are labelled, although 
the text is not visible from some viewing angles. 

The 2D place-holder figure shows two position-velocity 
maps, constructed by summing the data cube slice over the 
y- (left) and x- (right) axes. The ionized gas cloud occupies 
a very restricted region of velocity space, compared to the 
main body of the galaxy which has a clear rotation curve 


^ The full documentation for the Mayavi2 package is available at 



with a velocity amplitude of ~300 kms“^. The velocity of 
the gas cloud is very similar to that of the central AGN, 
consistent with the velocity fields shown in Figure [7o| The 
clear morphological and kinematic association between the 
gas cloud and the main body of the galaxy indicates that 
there is a strong link between these two physical components 
of J2240. 

The line-emitting gas in the cloud typically has line 
ratios and electron temperatures consistent with ionization 
due to AGN activity. The [GUI] peaks have large electron 
densities (ue ~125 cm“^) compared to the underlying gas 
in the cloud (ue < 50 cm“^); suggesting that gas may be 
compressed at shock fronts along the AGN ionization cone. 
The northern [O III] peak has the largest velocity disper¬ 
sion in the cloud (clear from the large range of gas velocities 
visible at this spatial location; see the left hand panel of 
the 2D figure and the X position - velocity view of the 3D 
figure). The high velocity dispersion combined with strong 
[Nelll] flux and large [NII]/Ha and [SII]/Hq: ratios sug¬ 
gests that this region is strongly shock excited. The two 
southern [O III] peaks have relatively small velocity dis¬ 
persions, large [OIII]/H/3 ratios and small [NII]/Ha and 
[SII]/Hq: ratios, indicating that the AGN is the dominant 
ionizing source in these regions. However, the detection of 
[Nelll] and Hell emission combined with electron temper¬ 
atures of ~20,000 K (especially in the peak closest to the 
central AGN) suggests that shock excitation is likely to con¬ 
tribute to the energy budget of these [O 111]-bright regions. 
The low velocity dispersions are consistent with shocks prop¬ 
agating mainly perpendicular to the line of sight, along the 
AGN ionization cone. 

The overlap region between the cloud and the main 
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body of the galaxy is characterised by a superposition of 
the disc and cloud rotation curves. Double peaked emission 
line profiles are observed but restricted to a small region 
in which the emission from both components is of approxi¬ 
mately equal intensity (see Figure]^. 

The main body of the galaxy has a compound emission 
signature which reveals both a narrow-line rotating compo¬ 
nent and a prominent broad component which is ubiquitous 
across the main body of the galaxy. Strong [Ne V] emission 
signposts the location of the AGN at the centre of the galaxy 
disc. The narrow emission component has similar line ra¬ 
tios and electron temperature to the line-emitting gas in the 
cloud, indicating that both are ionized primarily by the cen¬ 
tral AGN. The majority of the underlying broad component 
is also likely to be energized by the central AGN, whilst lo¬ 
calised regions of fast shocks can be identified by their large 
[NII]/Ha and [S II]/Ha ratios, strong [Ne III] and He H flux 
and large velocity dispersions. 


The elongation on the north eastern side of the galaxy 
is seen in both continuum and line emission, indicating that 
J2240 may have been morphologically disrupted during past 
interactions with a companion galaxy. Strongly redshifted, 
turbulent gas suggests the presence of outflowing material. 
Large [OHI]/H^, [NII]/Ha and [SH]/Ha ratios and an 
electron temperature of ~25,000 K provide additional ev¬ 
idence for significant shock excitation in the region. 


The [NelH] peak on north eastern side of the galaxy 
disc (see Section [3^ is a particularly interesting feature of 
J2240. It has a much larger [O HI] flux than the neighbour¬ 
ing regions of the galaxy, and its Hell flux is comparable 
to that of the central AGN. It has a peak velocity disper¬ 
sion of 500 kms“^, an electron temperature of ~ 20,000 K, 
and similar line ratios to the AGN. The strong warp in the 
galaxy’s [O III] emission profile appears to be driven pri¬ 
marily by strong emission from this concentrated region 
of turbulent, hot, dense gas. This hot spot may be en¬ 
ergized by an AGN-driven jet or outflow interacting with 
the interstellar medium. J2240 has a 1.4 GHz luminosity of 
1.08 X lO^^WHz”^, corresponding to a radio power of 6.8 


X10^® ergs ^ and a jet power of ' 


10^^ ergs ^ (using the 


Pjet — P 1.4GHz scaling relation from Gavagnolo et al. polol l. 
The localized nature of the shocked spot suggests that the 
interacting material is highly collimated, and therefore the 
hotspot is unlikely to be energised by an outflow alone. High 
resolution Very Long Baseline Interferometry (VLBI) obser¬ 
vations are required to resolve and further analyse the radio 
structure of this target. 


Although each of the features we have discussed have 
distinct physical properties and are clearly resolved spec¬ 
trally, the minimum separations between the features are 
1-2 resolution elements (where one resolution element is the 
FWHM of the PSF) and therefore they are not spatially re¬ 
solved. This does not impact our conclusions regarding the 
ionization sources and kinematics of J2240. Our 3D spectro¬ 
scopic analysis has indicated that decomposing the spectra 
into different kinematic components allows us to spectrally 
resolve distinct ionizing sources, even when they are super¬ 
imposed along the line of sight and therefore cannot be spa¬ 
tially resolved. 


5.2 Fraction of [O III] luminosity due to AGN 
activity 

The large [O HI] luminosity of J2240 makes it a very unusual 
object. The [O III] luminosity is both an order of magnitude 
larger than expected from its 22/im luminosity and an or¬ 
der of magnitude larger than that of typical Type-2 quasars 
at similar redshifts. However, our spectroscopic analysis has 
indicated that shock excitation is a non-negligible ioniza¬ 
tion source which can account for some of this discrepancy. 
We estimate the relative contribution of shock excitation to 
the [O HI] luminosity of J2240 by assuming that shock ex¬ 
citation is responsible for all component 2 emission in the 
cloud region of the galaxy and all component 3 emission to 
the north east of the AGN. This is a conservative estimate 
given that the AGN narrow line region extends across the 
entire [O 111]-emitting region of the system. We find that the 
AGN is responsible for at least 82 per cent of the galaxy’s to¬ 
tal [O HI] luminosity (this number rises to 86 per cent when 
correcting for extinction using an unreddened Balmer decre¬ 
ment of 3.1). The contribution from shocks is clearly non- 
negligible and reduces the discrepancy between the [O III] 
and 22/im luminosities by a factor of ~1.2, but is not large 
enough to account for the full order of magnitude discrep¬ 
ancy. This confirms that J2240 is a very unusual object, 
illuminated by the radiation field of an AGN which appears 
to have drastically reduced its accretion rate over less than 
the light crossing time from the central engine to the narrow 
line region. 


5.3 Origin of the Extended Emission Line Region 
around J2240: remnant of a quasar-driven 
outflow? 

To date, the majority of EELRs have been observed around 
three main classes of galaxies: the most powerful radio loud 
quasars at intermediate redshifts (0.15 < z <0.4; Stockton 


MacKenty|1987 Eu V Stockton|200^ , low redshift galax¬ 


ies with weak central AGN (z < 0.1; Keel et al. 2012a 
2014) and radio quiet luminous obscured quasars at z < 0.8 


^ Greene et al.|20TT Liu, Zakamska V Greene|2014 Hainline 


et al. 2014). The EELRs around the first two classes of galax¬ 


ies extend to tens of kpc and are morphologically distinct 
from their host galaxies, whereas ionized gas around radio 
quiet luminous obscured quasars is ubiquitous to < 10 kpc 
but rarely extends beyond this radius, and appears to be 
morphologically associated with the host galaxy. This di¬ 
chotomy suggests that the largest EELRs are either tidal in 
origin or outflow remnants ( Eu V Stockton|2009 Keel et al 
2014 ), whilst smaller EELRs are limited by the availability of 
gas to ionize within the host galaxy (Greene et al.|20TT Liu, 


Zakamska V Greene|2014 Hainline et al.|20T4 ). The EELR 


of J2240 extends over at least 22 x 35 kpc (Schirmer et al. 


2013) and is morphologically distinct from the host galaxy. 


and therefore clearly falls into the first group of EELRs. 

The EELRs around intermediate redshift radio-loud 
quasars are thought to be remnants of quasar-driven su¬ 
perwinds triggered by galaxy interactions (Eu & Stockton 


2009). These regions are ionized by current strong AGN ac¬ 


tivity. The EELRs are morphologically distinct from their 
central galaxies and many contain bright, dense clumps of 
gas illuminated by shocks propagating through the gas. The 
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gas clouds have high velocities but low velocity dispersions, 
and appear to be kinematically ordered on small scales but 
disordered on large scales ( Fu Stockton|2Q09 ). In contrast, 
EELRs around low redshift galaxies with weak central AGN 
are thought to be tidal debris illuminated by AGN ioniza¬ 
tion echoes. The clouds are often diffuse and filamentary, 
resembling tidal tails or radial shell/ring structures. The gas 
clouds have line ratios consistent with AGN ionization, but 
otherwise appear to be kinematically very stable, following 
the ordered rotation curves of the galaxies they surround 


(Keel et al. 2014). The objects in the Keel et al. (2014) 


sample were selected to have at least a factor of 3 discrep¬ 
ancy between the observed AGN power and the AGN power 
required to illuminate the gas in the galaxy, and therefore 
provide an interesting comparison sample for GBs. 

The morphology, ionization properties and kinematics 
of J2240 indicate that it is either a different type of object 
or at a different evolutionary stage to the objects in both 


the Keel et al. (2014) and Eu & Stockton (2009) samples. 


J2240 is neither radio loud nor a quasar host galaxy, but 
drives a radio jet with a synchrotron power of ~ 10^^ ergs“^ 
which interacts with the interstellar medium of the galaxy 
and produces a bright, dense, turbulent hotspot. The [O III] 
luminosity of the EELR exceeds that of the most luminous 
quasars at the same redshift. This indicates that the quasar 
phase of J2240 must have ended more recently than the 
light crossing time from the central en gine to the EELR 
10^ — 10^ years; Schirmer et al.|20T3 ). We find clear evi¬ 


dence for hot, dense, luminous knots of gas along the ioniza¬ 
tion cone of the AGN. There is clear evidence for morpholog¬ 
ical disturbance in both the stellar continuum and ionized 
gas emission of J2240, consistent with the hypothesis that 
merger activity is a defining feature of systems with promi¬ 
nent EELRs. The entire galaxy appears to be embedded in 
a very low surface brightness asymmetric halo (see Eigure 8 
of Schirmer et al.||2013 ), but the morphology of the ionized 
gas cloud suggests that it is unlikely to be composed of tidal 
debris. The morphology and [O III] luminosity of the EELR 
therefore suggest that J2240 may be a direct descendent of 
the luminous quasars observed in the |Eu &: Stockton (2009) 
sample. 

The kinematic properties of J2240, however, point to¬ 
wards a more complicated evolutionary sequence, as noted 
by Schirmer et al. (2013). The EELR gas clouds appear to 


follow the large-scale rotation of J2240, making their motion 
globally ordered. The velocity structure of the line profiles 
in the south western region of the galaxy suggest that the 
EELR gas clouds have low velocities with high velocity dis¬ 
persions. However, our multiple component analysis has in¬ 
dicated that the average velocity dispersion of the EELR gas 
is around 100 kms“^, and that the higher velocity dispersion 


components identified by Schirmer et al. (2013) originate 


from gas associated with the galaxy disc which is superim¬ 
posed in front of the cloud along the line of sight. We do not 
observe any high velocity gas in the cloud, perhaps because 
the gas in the cloud has very little motion in the radial direc¬ 
tion or because the high velocity components are very faint 
and therefore undetected. The overall kinematic structure 
of the EELR in J2240 appears to be approximately interme¬ 


diate between that of the Eu & Stockton (2009) and Keel 


EELR is the remnant of a quasar-driven outflow rather than 
being tidal in origin. One possible evolutionary scenario is as 
follows. Strong gravitational interaction with a companion 
galaxy funnelled large amounts of cold gas towards the cen¬ 
tre of the system, triggering a powerful quasar phase and si¬ 
multaneously launching a radio jet and a quasar-driven out¬ 
flow. A collimated outflow would naturally produce outflow¬ 
ing material with a coherent velocity structure and a strong 
morphological link to the central galaxy. The [O III]-bright 
knots along the AGN ionization cone are consistent with 
compression of outflowing material along fast shock fronts. 
The ordered rotation of the gas in the cloud is consistent 
with its close proximity to the central galaxy. A sustained 
quasar driven outflow would have the ability to remove the 
majority of the cold gas from the galaxy, starving the cen¬ 
tral engine and perhaps contributing to the rapid decrease 
in the accretion rate of the AGN. Molecular gas observa¬ 
tions are required to further constrain the mass, kinematics 
and spatial distribution of molecular gas across the entire 
system. 

It is worth noting that J2240 is significantly different 
from Hanny’s Voorwerp. The ionization parameter in the 
gas surrounding the central engine of IG 2497 is much lower 
than in J2240 (logU ~ —3.2 compared to \ogU > —2.6). 
The highest measured velocities in the Voorwerp are on the 
order of ~ 90kms“^ compared to J2240 where we mea¬ 
sure velocities up to 400 kms“^ in the EELR and up to 600 
kms“^ in the main body of the galaxy. The average electron 
temperature in the Voorwerp is ~ 13,500 K but tempera¬ 
tures as large as 25,000 K are measured in J2240. We also 
detect strong Hell emission in many [O 111]-bright clumps 
offset from the centre of J2240, which is not observed in 
the Voorwerp. All of these differences indicate that shock 
excitation is much more prominent in J2240 than in the 
Voorwerp - IG 2497 system. However, both systems show 
evidence for the onset of outflows corresponding closely to 
the rapid decrease in brightness of their central engines; sug¬ 
gesting that the evolutionary scenarios of the two systems 
may not be significantly different. 

Whilst the EELR of J2240 appears to be the remnant 
of a quasar-driven outflow, the GB selection criteria do not 
guarantee the presence of gas >10 kpc from the nucleus of 
the host galaxy. The angular extents of the 29 GBs presented 


by Schirmer et al. (2013) imply physical Petrosian radii of 
4.5 - 15 kpc, and 71% of the sample have Petrosian radii 
< 10 kpc, suggesting that gas reservoirs extending > 10 kpc 
may not be a defining characteristic of the GB phenomenon. 
Eollow up observations of a larger sample of GBs are re¬ 
quired to determine whether or not their EELRs have a 
common origin, and whether differences in the physical prop¬ 
erties of EELRs are driven primarily by intrinsic properties 
of the host galaxy or by the mechanism which triggers the 
rapid decrease in accretion rate of the central engine. 


5.4 Space density of quasar ionization echoes 

The space density of GBs (4.4 Gpc“^) is very low, suggest¬ 
ing that the progenitors of GBs must be extremely rare. We 
use the redshift-dependent quasar luminosity function de¬ 


et al. (2014) samples. 

Our observations of J2240 strongly suggest that its 


rived by Ross et al. (2013) to investigate the expected space 
density of galaxies representing a ~ 10^ year phase in the 
lifetimes of luminous quasars. We convert the [O III] lumi- 
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nosities of the GBs into approximate rest-frame z-band mag¬ 


nitudes using the conversion given by Reyes et ah (2008). 


All GBs in our sample have Mi < — 25. We integrate the 
quasar luminosity function from Mi = -25 to Mi = -oo at 
both z = 0.19 and z = 0.35, the boundaries of the redshift 
window for the lower redshift GB sample. We find that the 
expected space densities of luminous quasars at z = 0.19 and 
z = 0.35 are ~ 1300 Gpc“^ and ~ 3300 Gpc“^ respectively. 
We then use a simple timescale argument to determine the 
expected ratio of quasars to GBs. We assume that the life¬ 
times of luminous quasars are approximately equivalent to 
the Salpeter time ( the e-folding t ime for Eddington-limited 
black hole growth; Salpeter||l964 ) which is 4.2 x 10^ years 


(see e.g. Hopkins h Hernquist||2009 ). We also assume that 


GBs is 10^ years (e.g. 

Schirmer et al.|2013 Keel et al.|2014 

2012a Lintott et al. 

2009). We then multiply the derived 


quasar space densities by the ratio of the GB lifetime to the 
quasar lifetime to derive estimates for the expected space 
densities of GBs. We find expected space densities of 3.1 
Gpc“^ and 7.9 Gpc“^ at z = 0.19 and z = 0.35 respectively. 
The observed space density of GBs lies in between these two 
numbers and is therefore consistent with a scenario in which 
GBs are a common but short-lived stage in the evolution of 
the most luminous quasars {Mi < -25 at 0.3 < z < 0.7). 


6 SUMMARY AND CONCLUSIONS 

We have used Gemini-GMOS IFU data to conduct a detailed 
analysis of the 3D ionization state and kinematics of the 
GB galaxy J2240. The majority of the optical emission from 
this system is concentrated in strong, high equivalent width 
nebular emission lines (e.g. [GUI] A5007, [O III] A4959, H/3, 
Ha, [N H], [O H], [S H]) which overpower the underlying stel¬ 
lar continuum emission. The emission line profiles are very 
complex, indicating that several kinematic components con¬ 
tribute significantly to the observed emission signature of 
the galaxy. We use the parametric emission-line fitting code 
LZIFU (Ho et ah, in prep) to model the spectrum of each 
spaxel as a linear superposition of up to three Gaussian com¬ 
ponents with varying velocities and velocity dispersions and 
thus reconstruct the emission spectrum of each kinematic 
component individually. Our analysis of the [O HI] morphol¬ 
ogy, velocity field, ionization state, and electron temperature 
and density of each of the kinematic components has indi¬ 
cated that each of the components have distinct physical 
properties: 

• Gomponent 1 has the smallest average velocity disper¬ 
sion (cr ~ 100kms“^) and traces the bulk motion of ion¬ 
ized gas in the cloud. The majority of the cloud material 
lies at the systemic velocity of the galaxy and has a strong 
morphological and kinematic link with the galaxy’s cen¬ 
tral engine. The remainder of the cloud material appears 
to be co-rotating with the main body of the galaxy but 
with a shallower rotation curve than that of the disc ma¬ 
terial. The strongest [O HI] emission comes from clumps of 
gas which lie along the ionization cone of the central AGN. 
The [OHI]/H/3, [NII]/Ha and [SH]/Ha ratios of the gas in 
the cloud indicate that the AGN is the dominant ionizing 
source. The average electron temperature is Te ~14,700 K, 


consistent with previously measured values for AGN narrow 
line regions. 

• Gomponent 2 has a slightly larger average velocity dis¬ 
persion (a ~ 150kms“^) and traces the bulk motion of gas 
associated with the main body of the galaxy, as well as 
shocked material in the cloud. The large scale gas motions 
reveal a rotation curve with an amplitude of z; ~ 300 km s“^ 
across the main body of the galaxy. The centre of this rota¬ 
tion curve is coincident with the peak of the galaxy’s con¬ 
tinuum emission which is also a [Ne V] source and is there¬ 
fore likely to be the central AGN. The diagnostic line ra¬ 
tios indicate that the AGN is again the dominant ionizing 
source. However, [NelH] and Hell detections in dense, tur¬ 
bulent clumps of gas in the cloud indicate that shock exci¬ 
tation also contributes significantly to the ionization signa¬ 
ture of this component. The average electron temperature 
is Te -12,900 K. 

• Gomponent 3 traces very turbulent gas 
(a — 300 — 800kms“^) in the main body of the galaxy, 
ionized by fast shocks and the hard ionizing radiation field 
of the central AGN. The average electron temperature in 
this component is Te — 19,300 K. There is some evidence 
for ordered rotation, but there are many pockets of gas 
whose velocities deviate significantly from the underlying 
rotation curve. The diagnostic line ratios indicate that 
the majority of this component is likely to be excited by 
the hard ionizing radiation field of the AGN. However, 
[NelH] and Hell detections in regions with large [NII]/Ha 
and [SII]/Ha ratios and high electron temperatures 
(Te — 20,000 K) again indicate that fast shocks are a 
prominent ionizing source in this component. We identify 
a turbulent hotspot to the north east of the AGN whose 
properties are consistent with a radio-jet shock-heating the 
interstellar medium. 


Our analysis confirms that J2240 is a very unusual ob¬ 
ject, with disturbed morphology, complex kinematics and 
several ionizing sources. Our results are both qualitatively 
and quantitatively consistent with the findings of |Schirmer| 


et al. (2013), but provide a much more detailed insight into 


the principal ionizing sources of J2240, allowing us to esti¬ 
mate the contribution of the AGN to the total [O III] lu¬ 
minosity. We find that the AGN must be responsible for at 
least 82 per cent of the total [OIII] luminosity of J2240. 
Shock excitation is a non-negligible ionizing source but is 
not significant enough to resolve the discrepancy between 


the [O HI] and 22/im luminosities of J2240 (see Schirmer 
et al.|[^13| . This confirms that the EELR in J2240 is 


quasar ionization echo. 

The evolutionary history of J2240 remains a puzzle. Al¬ 
though the EELR of J2240 has slightly different properties 
to similar regions observed around radio-loud quasars, our 
observations indicate that the cloud may be the remnant of 
a quasar-driven outflow, triggered by previous a galaxy in¬ 
teraction and now illuminated by a quasar ionization echo. 
VLBI and molecular gas observations are required to con¬ 
strain the morphology of the radio emission and the mass, 
kinematics and spatial distribution of molecular gas in the 
system. 3D spectroscopic analysis and deep imaging of a 
larger sample of GBs will reveal whether these objects share 
similar properties to J2240, and ultimately will allow us to 
determine the evolutionary history of the GB phenomenon. 

































22 R. L. Davies, M. Schirmer, J. E. H. Turner 


7 ACKNOWLEDGEMENTS 

The authors would like to thank the referee for their helpful 
comments which improved the clarity of this manuscript. RD 
acknowledges the support of an Australian Gemini Under¬ 
graduate Summer Studentship from the Australian Astro¬ 
nomical Observatory funded by Astronomy Australia Ltd. 
RD would like to thank Geoff Bicknell for valuable discus¬ 
sions regarding radio jets and AGN-driven outflows. Based 
on observations obtained at the Gemini Observatory (ac¬ 
quired through the Gemini Science Archive and processed 
using the Gemini IRAF package), which is operated by the 
Association of Universities for Research in Astronomy, Inc., 
under a cooperative agreement with the NSF on behalf of 
the Gemini partnership: the National Science Foundation 
(United States), the National Research Gouncil (Ganada), 
GONIGYT (Ghile), the Australian Research Gouncil (Aus¬ 
tralia) , Ministerio da Giencia, Tecnologia e Inovacao (Brazil) 
and Ministerio de Giencia, Tecnologfa e Innovacion Produc- 
tiva (Argentina). Ureka is a collaboration of Space Telescope 
Science Institute and Gemini Observatory. We thank Bryan 
Miller for his contributions to the ifudrgmos package used 
for processing our data. 


REFERENCES 

Alexander D. M., Swinbank A. M., Small I., McDermid R., 
Nesvadba N. P. H., 2010, MNRAS, 402, 2211 
Allen M. G., Groves B. A., Dopita M. A., Sutherland R. S., 
Kewley L. J., 2008, ApJS, 178, 20 
Allington-Smith J. et ah, 2002, PASP, 114, 892 
Asmus D., Gandhi P., Smette A., Honig S. F., Duschl W. J., 
2011, A&A, 536, A36 

Baldwin J. A., Phillips M. M., Terlevich R., 1981, PASP, 
93, 5 

Bertelli G., Bressan A., Ghiosi G., Fagotto F., Nasi E., 
1994, A&AS, 106, 275 
Bessell M. S., 1999, PASP, 111, 1426 
Gappellari M., Gopin Y., 2003, MNRAS, 342, 345 
Gappellari M., Emsellem E., 2004, PASP, 116, 138 
Gavagnolo K. W., McNamara B. R., Nulsen P. E. J., Garilli 
G. L., Jones G., Birzan L., 2010, ApJ, 720, 1066 
Dopita M. A., Sutherland R. S., 1995, ApJ, 455, 468 
Evans I., Koratkar A., Allen M., Dopita M., Tsvetanov Z., 
1999, ApJ, 521, 531 

Eischera J., Dopita M., 2005, ApJ, 619, 340 
Eu H., Stockton A., 2009, ApJ, 690, 953 
Gaskell G. M., Eerland G. J., 1984, PASP, 96, 393 
Girardi L., Bertelli G., Bressan A., Ghiosi G., Groenewegen 
M. A. T., Marigo P., Salasnich B., Weiss A., 2002, A&A, 
391, 195 

Girardi L., Bressan A., Bertelli G., Ghiosi G., 2000, A&AS, 
141, 371 

Gonzalez Delgado R. M., Gerviho M., Martins L. P., Lei- 
therer G., Hauschildt P. H., 2005, MNRAS, 357, 945 
Goudfrooij P., Bohlin R. G., Mafz-Apellaniz J., Kimble 
R. A., 2006, PASP, 118, 1455 
Greene J. E., Zakamska N. L., Ho L. G., Barth A. J., 2011, 
ApJ, 732, 9 

Greisen E. W., Galabretta M. R., Valdes E. G., Allen S. L., 
2006, AVA, 446, 747 


Groves B. A., Dopita M. A., Sutherland R. S., 2004, ApJS, 
153, 9 

Hainline K. N., Hickox R., Greene J. E., Myers A. D., Za¬ 
kamska N. L., 2013, ApJ, 774, 145 
Hainline K. N., Hickox R. G., Greene J. E., Myers A. D., 
Zakamska N. L., Liu G., Liu X., 2014, ApJ, 787, 65 
Hamuy M., Suntzeff N. B., Heathcote S. R., Walker A. R., 
Gigoux R, Phillips M. M., 1994, PASP, 106, 566 
Ho I.-T. et ah, 2014, MNRAS, 444, 3894 
Hook 1. M., Jprgensen L, Allington-Smith J. R., Davies 
R. L., Metcalfe N., Murowinski R. G., Grampton D., 2004, 
PASP, 116, 425 

Hopkins P. E., Hernquist L., 2009, ApJ, 698, 1550 
Hopkins P. E., Quataert E., 2010, MNRAS, 407, 1529 
Ichikawa K., Ueda Y., Terashima Y., Oyabu S., Gandhi P., 
Matsuta K., Nakagawa T., 2012, ApJ, 754, 45 
Keel W. G. et ah, 2012a, MNRAS, 420, 878 
Keel W. G. et ah, 2012b, AJ, 144, 66 
Keel W. G. et ah, 2014, ArXiv e-prints 
Kewley L. J., Dopita M. A., 2002, ApJS, 142, 35 
Kewley L. J., Dopita M. A., Leitherer G., Dave R., Yuan 
T., Allen M., Groves B., Sutherland R., 2013, ApJ, 774, 
100 

Kewley L. J., Dopita M. A., Sutherland R. S., Heisler G. A., 
Trevena J., 2001, ApJ, 556, 121 
Khan E. M., Berentzen L, Berczik P., Just A., Mayer L., 
Nitadori K., Gallegari S., 2012, ApJ, 756, 30 
Levesque E. M., Richardson M. L. A., 2014, ApJ, 780, 100 
Lintott G. J. et ah, 2009, MNRAS, 399, 129 
Liu G., Zakamska N. L., Greene J. E., 2014, MNRAS, 442, 
1303 

Liu G., Zakamska N. L., Greene J. E., Nesvadba N. P. H., 
Liu X., 2013, MNRAS, 430, 2327 
Markwardt G. B., 2009, in Astronomical Society of the Pa¬ 
cific Gonference Series, Vol. 411, Astronomical Data Anal¬ 
ysis Software and Systems XVHI, Bohlender D. A., Du¬ 
rand D., Dowler P., eds., p. 251 
Massey P., Strobel K., Barnes J. V., Anderson E., 1988, 
ApJ, 328, 315 

Mullaney J. R., Alexander D. M., Eine S., Goulding A. D., 
Harrison G. M., Hickox R. G., 2013, MNRAS, 433, 622 
Nicholls D. G., Dopita M. A., Sutherland R. S., 2012, ApJ, 
752, 148 

Nicholls D. G., Dopita M. A., Sutherland R. S., Kewley 

L. J., Palay E., 2013, ApJS, 207, 21 

Osterbrock D. E., Eerland G. J., 2006, Astrophysics of 
gaseous nebulae and active galactic nuclei 
Osterbrock D. E., Eulbright J. P., Martel A. R., Keane 

M. J., Trager S. G., Basri G., 1996, PASP, 108, 277 
Perez-Montero E., Hagele G. E., Gontini T., Dfaz A. L, 

2007, MNRAS, 381, 125 

Poppett G. L., Allington-Smith J. R., 2010, MNRAS, 404, 
1349 

Ramachandran P., Varoquaux G., 2011, GSE, 13, 40 
Reyes R. et ah, 2008, AJ, 136, 2373 

Rich J. A., Dopita M. A., Kewley L. J., Rupke D. S. N., 
2010, ApJ, 721, 505 

Rosario D. J., Shields G. A., Taylor G. B., Salviander S., 
Smith K. L., 2010, ApJ, 716, 131 
Ross N. P. et ah, 2013, ApJ, 773, 14 
Salpeter E. E., 1964, ApJ, 140, 796 
Schawinski K. et ah, 2010, ApJ, 724, L30 



J224024-1-092748: Quasar ionization echo 23 


Schirmer M., Diaz R., Holhjem K., Levenson N. A., Winge 
C., 2013, ApJ, 763, 60 

Stockton A., Fu H., Canalizo G., 2006, New A Rev., 50, 
694 

Stockton A., MacKenty J. W., 1987, ApJ, 316, 584 
Tanaka T., Haiman Z., Menou K., 2010, AJ, 140, 642 
Tody D., 1986, in Society of Photo-Optical Instrumentation 
Engineers (SPIE) Conference Series, Vol. 627, Instrumen¬ 
tation in astronomy VI, Crawford D. L., ed., p. 733 
Tody D., 1993, in Astronomical Society of the Pacific Con¬ 
ference Series, Vol. 52, Astronomical Data Analysis Soft¬ 
ware and Systems II, Hanisch R. J., Brissenden R. J. V., 
Barnes J., eds., p. 173 

Veilleux S., Osterbrock D. E., 1987, ApJS, 63, 295 
Vogt E. P. A., Dopita M. A., Kewley L. J., 2013, ApJ, 768, 
151 


APPENDIX A: MODIFICATIONS TO CMOS 
IRAF PACKAGE 

We reduced our data using a modified version of the Gemini 
CMOS IRAF package, as described briefly in Section |2.2| 
Here we include a detailed description of our modifications. 


A1 Bias and overscan 

To avoid errors due to charge-transfer contamination in the 
overscan region, the bias zero-level was determined with the 
gfreduce.biasrows parameter set to “3:64” for GMOS-S and 
“6:46” for GMOS-N with the EEV deep-depletion (EEV- 
DD) CCDs, selecting a small, unilluminated section at one 
end of the overscan region. 


A2 Bad columns 

For our GMOS-S data, we used gemarith to create an inte¬ 
ger bad pixel mask from a science exposure and imreplace 
to flag pixels that correspond to obvious bad columns (and 
a small, associated glow region), based on inspection. Sep¬ 
arate masks were produced for science, standard star and 
calibration exposures, as the presence or extent of some bad 
columns was observed to vary between exposures of differ¬ 
ent length and brightness. The masked pixel values were re¬ 
placed with a low-order (5th-32nd, depending on the image 
structure) fit to each row (spectrum), using a new wrapper 
task, gemfix, which has options to use IRAF’s fitld or fixpix. 

For the GMOS-N EEV-DD detectors, no bad columns 
were apparent, so no bad pixel masks were applied prior to 
cosmic ray rejection. 


A3 Scattered light 

Scattered light in the background of flat-field and scientific 
spectra was modelled by fitting a 5th- to 9th-order surface 
to the values of unilluminated image rows between blocks of 
fibres, over each CCD amplifier. This model was then sub¬ 
tracted from the data, providing an order-of-magnitude im¬ 
provement in the nominally-unilluminated zero-levels. The 
contamination in question is present at the level of up to ~ 15 
per cent over the useful range of the data (deteriorating to 


~100 per cent at the faint extremes). This processing was 
achieved using the ifudrgmos scripts fndblocks and gfbkgsub. 
Flat-field spectra were extracted both before and after this 
procedure, first to determine the fibre and gap locations, as 
input to fndblocks, and then with the scattered light removed 
for more accurate results. 


A4 Wavelength calibration 

The wavelength calibrations for our blue and green chan¬ 
nels are derived from the day-time CuAr spectra provided 
by Gemini. Only the ~20-40 brightest lines (with >100 
counts) were used, minimizing any centroid errors due to 
charge trailing from faint pixels during readout, which we 


have observed in previous CMOS datasets (see Goudfrooij 
et al.|2006| |. 


For the green channel, each exposure’s wavelength zero 
point was adjusted to account for possible flexure with re¬ 
spect to the arc, by measuring 13 lines from |Osterbrock| 
et al. (1996) in the sky spectrum using rvidlines and editing 


the FITS CRVALl keywords accordingly; this resulted in 
an RMS absolute correction of 0.5A and an RMS correction 
relative to the first exposure of 0.4A. 

For the blue channel, the only sky line detected ([N I] 
at 5199 A) had too low a S/N in two exposures to provide 
accurate flexure correction, so none was applied, but our 
low-S/N measurements indicate an RMS difference between 
exposures of 0.3 A (<0.lA between those with a strong de¬ 
tection), consistent with the moderate range of elevation 
(55-60°) and rotation (45°) of the observation. 

For the red channel, the provided CuAr spectrum is very 
sparse, so 17 strong sky lines were instead measured directly 
in each object exposure and compared with |Osterbrock et al.| 
(1996) to derive wavelength solutions. These solutions are 
consistent between fibres, unaffected by the much weaker 
object flux. Only two sky lines coincide with (still several 
times fainter) line emission from the source. Some ghost lines 
at the blue end of the red slit were ignored and do not overlap 
with our range of scientific interest. 


A5 Flat fielding 

A normalized flat field spectrum was produced using a ver¬ 
sion of gfextract that we have re-written previously to ac¬ 
count for wavelength distortions (the original version divides 
spectra with slightly differing wavelength solutions by a sin¬ 
gle vector of fitted continuum values). This change avoids 
minor artificial roll-off in flux towards the corners of the 
stack of extracted spectra and, in particular, eliminates ~5 
per cent discontinuities between the two halves of the IFU 
field. When fitting the continuum, our algorithm produces 
a steep drop at each extreme of the spectrum, so we use the 
gf response, sample parameter to restrict the fit to the useful 
wavelength range, plus the first and last pixels to constrain 
the ends. 

Our flats were generated using only GCal quartz- 
halogen exposures, omitting the optional twilight exposures 
(gfresponse.skyimage) that are often used to account for any 
slow turn-down in GCal intensity over the field. It is well 
known that gravitational flexure can lead to variations in 


the throughput of optical fibres (eg. Poppett & Allington- 
Smithll^TO ), depending on their individual stresses. When 
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processing GMOS-IFU data using twilight exposures, which 
are taken at a different pointing from the target, we have 
encountered some number of artihcahy bright or dark spec¬ 
tra, consistent with hexure effects, which we do not see using 
only contemporaneous GCal hats. Moreover, the 5x7 arcsec 
IFU is only a fraction of the size of the 5.5 arcmin GMOS 
imaging held, over which GGal hatness is much more of an 
issue, and we do not observe measurable roll-off over the IFU 
held. Even when test-processing twilight spectra in place of 
scientihc exposures, the spatial residuals do not exhibit a 
low-order trend at the detectable level of a few percent. 

A6 Cosmic ray cleaning 

Gosmic rays were identihed using the spectroscopic version 
of L.A. Gosmic for IRAF (via the gscrspec wrapper script 
in ifudrgmos), with a modihcation to the imcalc step at line 
277 to avoid artifacts in the presence of negative values. We 
then re-cleaned the applicable pixels using our new gemfix 
script with its fixpix option, providing slightly more accu¬ 
rate replacement values than L.A. Gosmic’s running median. 
Pixels hagged as cosmic rays are also later excluded during 
spatial interpolation by gfcuhe. 

A7 Flux calibration 

The standards taken for this programme in photometric con¬ 
ditions were used for absolute spectrophotometric hux cali¬ 
bration. These were processed similarly to science exposures 
but with only light (~ lOcr) cosmic ray cleaning, which is suf- 
hcient to avoid a bias due to cosmic rays while avoiding spu¬ 
rious rejection caused by L.A. Gosmic not being optimized 
for the strong modulation in bright hbre IFU spectra. The 
gap between the two brightest blocks of hbres was omitted 
when htting scattered light for the blue channel, due to the 
strength of the PSF wings in that nominahy-unihuminated 
region. For the red channel, the sky lines in the longer scien¬ 
tihc exposure were used for wavelength calibration. A simple 
hux summation over the IFU object hbres (with gfapsum) 
was used to extract the standard spectra. For the green chan¬ 
nel, a non-photometric standard taken close in time to the 
science data was used to derive an initial response function, 
whose absolute level was then corrected using a standard ob¬ 
tained in photometric conditions 2 months later (both giv¬ 
ing very close results). We conhrm from our standard PSFs 
in reconstructed images that a negligible proportion of hux 
was lost outside the area of the IFU held, so we have not 
corrected for this. 

Our hnal spectra were compared with tabulated 
huxes provided in the gmos$calih (for EG 131) or oned- 
stds$spec50cal IRAE directories, adding some interpolated 
bandpasses (with the hux in the surrounding bins varying 
by ^ 5 per cent) where necessary to constrain the ends of 
the ht, mainly for the red channel. After applying the re¬ 
sulting sensitivity function to the science data, we used a 
new gfcube option to convert the hux per hbre to units of 
erg cm“^ s~^ arcsec”^ in our data cubes. 

No compensation was applied for slightly different vari¬ 
ations in quantum efficiency with wavelength between the 3 
individual GMOS GGDs, as no measurements are available 
to date for the GMOS-N EEV-DD detectors and the pro¬ 
cessing functionality is still under development by Gemini. 


Eor the GMOS-S GGDs we expect resulting discontinuities 
in hux at the level of ~2 per cent at the chip gap wave¬ 
lengths (B. W. Miller, private communication), which are 
not obvious in the data and would not signihcantly alter 
our results. 


A8 Reconstruction and combination of data cubes 


Our calibrated spectra were resampled to 3-dimensional 
data cubes with gfcube, using its built-in compensation 
for atmospheric dispersion (modelled using the SLALIB li¬ 
brary), which we added previously for a different project. 
We have also added an option to propagate approximate 
variance alongside the science array and used the resulting 
error estimates in our analysis. These errors are approximate 
because of the statistical correlation between adjacent pix¬ 
els generated by resampling lEU hbres onto a hner spatial 
grid in the data cubes. To avoid computationally-expensive 
propagation of covariance arrays, we simply upscale the vari¬ 
ance by a factor of \/T2 (the ratio of spatial pixel densities) 
such that noise estimates will be correct when averaging over 
more than ~ 3 pixels spatially. 

The data cubes corresponding to individual exposures 
were combined using our aforementioned PyFU Python 
scripts. We have added new options to facilitate regis¬ 
tration of nebulous structure via spatial cross-correlation 
{pyfalign.method=‘^correlate”) and to inspect the resulting 
alignment more easily {pyfmosaic.separatelyes), as well as 
propagating variance automatically. 

A new PyEU script, pyflogbin, was used to resample 
the co-added data cubes (and their variance) in bins of con¬ 
stant AA/A, for analysis, following the convention of Greisen 
et al.| ( [200^ , eq. 5; this is conveniently similar to the linear 
input WGS but differs from IRAE’s representation. This re¬ 
sulted in constant radial velocity increments of 28, 22 and 
11 kms“^ in the blue, green and red channels respectively. 








